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Penetration of Gamma Radiation From a Plane 
Monodirectional Oblique Source’ 


Martin J. Berger 


The penetration of gamma radiation from a plane monodirectional oblique source is 
solved by the moment method of Spencer and Fano. The novel feature of the treatment is 
the systematic investigation and exploitation of a function-fitting technique, which consists 
of approximating the desired radiation-flux function, f(7), by an expression of the form 
~A,o( Br , the parameters A, and B, being determined from the knowledge of the flux 
moments. 

A discussion is given of the criteria for the proper choice of the “weight function,”’ 
o(c), and of the accuracy of the fitting procedure. The results of a sample calculation are 
presented for an 0.66-Mev source in water. 


1. Introduction 


One of the most effective techniques for calculating the diffusion of gamma rays in infinite 
media is the moment method [1],? i. e., the determination of the spatial and angular moments 
of the radiation flux and subsequent calculation of the flux-distribution function by moment- 
fitting. This technique has been successfully applied to the calculation of radiation from plane 
and point isotropic sources and from monodirectional plane sources emitting radiation per- 
pendicularly to the source plane [2]. The purpose of this work is to extend the method to the 
determination of the propagation of radiation from monodirectional ob/ique plane sources, 1. e., 
sources emitting photons at a specified obliquity angle with respect to the source plane. In 
the problems treated previously, analytical asymptotic theory was available describing the 
radiation field at great distances from the source [3]; this facilitated the reconstruction of the 
flux funetion from its moments. For the problem of the oblique source, reliable guidance in 
the form of a corresponding analytical asymptotic theory is lacking, so that increased care must 
be given to the moment-fitting procedures. 

Although specific application of the method described in this paper will be made to the 
diffusion of gamma radiation, the technique used should also be useful for determining the 
propagation of other types of radiation, e. g., neutrons, from plane oblique sources. 

We consider an infinite homogeneous scattering medium that contains, in the plane z=0, a 
radiation source emitting 1 photon per square centimeter per second. The source photons are 
all emitted with energy /, and in directions inclined at an angle @ with respect to the normal 
to the source plane. (The distribution of the azimuthal directions with which the photons 
are emitted may be considered arbitrary, being irrelevant for the propagation as a function 
of distance from the source plane.) Our task is to determine the energy flux, fF“? (/,2)d&, 
at a distance 2 from the source plane, due to photons in the energy range (/,4+-dE) as a 
function of the obliquity @ The energy flux is defined as the flow of energy through an iso- 
tropic detector with unit cross-sectional area. The energy, /, will be expressed in million 
electron volts and the distance, 2, in units of the mean free path of the primary radiation. 

Unscattered radiation will be treated separately because its direct calculation is trivial. 
In order to facilitate the moment-fitting procedures, separate consideration will also be given to 
singly-scattered radiation. This separation is desirable because its flux distribution differs 
considerably from that of the multiply-seattered radiation. 


1 This work was supported by the Office of Naval Research and the Armed Forces Special Weapons Project 
? Figures in brackets indicate the literature references at the end of this paper 
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Section 2 discusses the moment-fitting procedures appropriate to the problem. In section 
3 specific application of these methods is made to radiation with a source energy of 0.66 Mey 
(cesium—137 radiation) in a scattering medium of water. Various energy spectra and spac: 
distributions are presented as functions of the obliquity angle. In section 4 the accuracy o|! 
the results is investigated by the comparison of various methods of moment-fitting and by 
comparisons with other related calculations and with an experiment. Details and subsidiary 
questions are given in appendixes. 

In appendix A, equations are derived for the calculation of the flux of singly-scattered 
radiation. In appendix B it is explained how the flux moments are calculated bv the Spencer- 
Fano method; in particular, it will be pointed out how the moments required for the treatment 
of the monodirectional oblique plane source radiation can be obtained from a certain set of 
standard moments computed on the NBS automatic computer (SEAC). Appendix C presents 


details of the moment-fitting procedure 
2. Moment-Fitting Methods 
2.1. Introductory Remarks 


We assume that the first eight spatial moments of the energy flux, 


Z"(E a) 2" Fo) (Ff z)dz, (1) 


are known. (It is assumed that F‘@” stands for the flux of singly-scattered or multiply-seattered 
radiation, whichever we may want to determine. For details concerning the computation of 
the moments, see appendix B. 

Although a well-behaved function is determined uniquely if all its moments are known 
the reconstruction of a function from only eight moments requires some finesse. It is desirable 
to incorporate prior knowledge concerning the flux that is available from related analytical 
calculations, physical intuition, or experience with similar problems, into a so-called weight 
function, which serves as the initial approximation. Successive improved approximations 
are then found by the use of the moments. 

In order to choose this weight function properly, we must give some thought to the ex- 
pected spatial dependence of the flux function. Assuming that OSa@<_w/2, 1. e., that the 
photons are released by the source in the direction of positive z, we anticipate that the flux of 
scattered radiation will be a unimodal asymmetric function of 2 with a peak at a value of 
somewhat larger than zero, and decreasing more rapidly on the negative than on the positive 
side of the peak. The asvmmetry of the funetion will depend on the obliquity, being a maxi- 
mum when a=0, and zero when a=7/2. Furthermore, it can be expected that at large dis- 
tances from the source plane the flux has a dominantly exponential behavior as function of 

To find a weight function embodying these features is difficult. The task is made easier 


if the flux is split into a symmetric component /°?(/,2) and an antisvmmetric component 
rs Ie.=), so that 

PO (E, f@ (E,z)+-9@(E,2 (2a) 

Po (E,—z)=f @(E,z)—g@ (E£,2), (2b) 

where 2 denotes |z!. In this manner we can account for the expected peak without incorporat- 


ing this feature into {6 and g®@, which may be chosen to be more or less smooth, monotonic 


asymptotically exponential functions. There is another advantage inherent in this splitting 
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chnique. We can fit four even moments to determine f‘® and four odd moments to deter- 
mine g@. This is much easier than fitting eight moments simultaneously. There is a slight 
disadvantage also; the flux cannot be determined very accurately far behind the source be- 
cause for large negative values of z the functions f‘” and g‘@ will have nearly the same value, 
ind their difference will be rather inaccurate. This is not a serious drawback because the 
small amount of radiation far behind the source usually is not of much interest. 


2.2. Polynomial Expansions 


We shall now discuss the method of polynomial expansion used by Spencer and Fano 
1] for the solution of the plane source with a=0 and explain why this approach cannot readily 
be carried over to the treatment of the source with a0. 

Spencer and Fano also made use of the splitting technique and set 


u 


: ‘ ¢ ; . — 9 
f(E,2=5- Doe (L)U (wz) (3a) 
— wze~™ 2. — ‘- 
g (E, 2) a — ay 4(L)V,(wz), (3b) 
~- n U 
where z denotes |2'; e~“?/2 and wze~“?/2 are the respective weight functions. 
: (—1)"/ 0 ae 2, (atj l 
U’,(2) : ( —1) ye gh me "ss (4a) 
2" O: 0 j (n—))! 
and 
9 O 2 
V’,(2) = ( - 1) U,(2) (4b) 
2(n+1)\ Oz 


(n=0,1,2,°° +) 


are orthogonal polynomial systems associated with these weight functions. The coefficients 
c,(f2) and d,(/2) are related to the moments as follows, the nth coefficient involving only 


moments of order not exceeding n. 





7 og “i { 1) “ 
c,(E)=> . mae (WZ) (oa) 
a  ) ai | 
, et fnt+l1\(-1! OS 
d(E)=>, (wz)?-! (5k 
HG Jeaj—ne" ' 


The validity of such a polynomial expansion depends on the correct choice of the exponent w 
in the weight functions. It can be shown by physical arguments and by related asymptotic 
penetration theory that a satisfactory choice for the plane source with obliquity zero is 


w= (u,/uo), Where wo is the linear narrow-beam absorption coefficient of the primary radiation 
and uw, that of the most penetrating component of the radiation. (This is usually, though not 
always, radiation with the primary energy.) With a correctly chosen parameter, the poly- 
nomial expansion will converge rapidly for wz up to ~15, so that only a few terms requiring 


the use of the first few moments have to be fed into the expansion. 

The parameter w is difficult to determine for the plane oblique source. Clearly, the radia- 
tion component with the highest energy (usually the component with the smallest attenuation 
coeflicient) will not govern the deep penetration because it must travel a distance approx- 
imately equal to 2/cos @ in order to get a distance 2 away from the source. On the other 
hand, radiation that has been turned parallel to the z-axis travels away from the source plane 
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on the shortest possible path but has lost energy, so that its attenuation is increased. On 
would expect that the deep penetration is controlled by the small amount of the radiatio: 
that has been turned approximately parallel to the z-axis in a number of collisions such tha 
it has lost very little energy. It has been found by numerical experimentation that it is diffi 
cult to obtain reliable polynomial expansions by adjusting the parameter w on the basis 
of convergence considerations. These difficulties can be overcome, however, by another 
approach [4]. 
2.3. Function Fitting 


If eight moments are available and an even-odd split is used, we set 


f@(E,2)=A\( Fb, a)¢o[2B,(E,a)|+ Al EL,a)o[2B.( F,a)| (Ga 
qe (E,z)=C (Ea ¥|2-D,\(E,a)| t Cha ¥(2D.(F,a)), (6b 
where $(z) and ¥(z) are suitable weight functions with the desired asymptotic and svmmetry 
properties: 
of } Ta | } laa 
¥ y } (7b 
yg re = }e , . (Na) 
7 _ >for large 2, 
¥(z)= Wz) (Sb 


where g(z) and ¥(z) are slowly varying functions of 2. 

The quantities A, and (; are intensity (normalization) parameters, and the 4, and D) 
are scale parameters. The introduction of location parameters is not necessary because the 
peak in the flux function will be determined by the superposition of svmmetrie and antisvm- 
metric solutions. The parameters 2; and J);, analogous to the parameter w in the method of 
polynomial expansion, will be effectively determined from the moments themselves, so that no 
“inspired guess”’ is necessary. 


We require that 


HK «= 37°C (E.a "+1 =D (EB.a)|d2. Ob 


The solution of these equations is described in appendix C. The important aspect of the solu- 
tion is the fact that there are restrictions on the class of allowed weight functions @(2) and W(< 
In order to obtain the desired asymptotic behavior, it is required that the scale parameters 
B, and D, are real-valued; otherwise, the weight functions would turn into oscillatory fune- 
tions, which is inconsistent with the expected smooth behavior of the radiation density function. 
As shown in appendix C, two pairs of dimensionless parameters may be introduced to 
characterize the relation between the function to be fitted and the proposed weight functions: 


ms M, M, A Mts MM 
om, m, M; , mm, AM? , 
1) 
mi MiM » mi MEM, 
Minn, M mym-, M 
where M,=2"(E.a),m ="o(z)dz, and Mo,44= hth 2) 2. 
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The weight function ¢(z) belongs to the allowed class (i. e., the scale parameters 4, are real) 
if, and only if, the point (s,f) lies in an “allowed region”’ of the s-t plane, indicated by the shaded 
area in figure | and specified in appendix C. ‘To determine whether ¥(2) is allowable, we replace 
the point (s,f) by (s,¢) and use the same diagram. 

It has been found in this investigation and in other related work, that considerable searching 
s often necessary to get an allowed weight function, but when there are several such functions 
they yield moment-fitting results that are in close agreement. Thus the function-fitting 
technique seems to operate in a sort of all-or-nothing fashion, and differs in this respect from 
the polynomial expansion method. 

The symmetric and asymmetric weight functions have been related as follows: 


¥(z)=z¢(2). (12) 


0, the antisymmetric component of the scattered flux vanishes and the 


This implies that at 
rather than some other antisymmetric 


total scattered flux is continuous.’ The choice of 2, 
function of 2 as the proportionality factor in eq (12) is somewhat arbitrary but recommends 
itself by its simplicity and is consistent with what is known about the analytical behavior ol the 


flux carried by photons with energies close to the source energy. 





———$—$—$— LT 
I unscatter ! »btained as space integral over the source function and the flux of the nth order of scattering as a space integral over 
flu the n—Ist order of scattering. The source function is a delta functior rhe unseattered flux thus has a simple discontinuity at the 
vhereas t flux of higher orders of scattering is continuou 
, 
1 
al 
>} 
toh 
6} 
F 
| 
4} 
OL - - ~ ~ ——— 
4 € 7 8 > 3 4 C 6 1.7 1.8 1.9 


Figure |. Diagram indicating whether a weight-function is “permissible” for function-filting. 


<}< but only a portion of t! revior shown. A weight-function 
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3. Numerical Application: Radiation from a 0.66-Mev Plane-Oblique Source 
in Water 


A sample calculation, based on the methods outlined in section, 2 was carried out for a 


source of 0.66-Mev photons (radiation from cesium-137) in water. 

The flux moments needed as input data were calculated on the SEAC at energies ranging 
from 0.66 to 0.134 Mev. At lower energies, spectra were obtained through an extrapolation 
procedure based on two characteristic features of gamma-ray diffusion: 

(1) At very low energies an equilibrium spectrum prevails in the sense that the flux can 
(at least, approximately) be factored into the product of a space-dependent function and an 
energy-dependent function. This has been shown to be the case for radiation from point and 
plane isotropic sources [2, 3], and there is no reason to doubt that it will also hold for plane 
oblique sources. 

(2) The shape of the energy spectrum at low energies is nearly independent of the source 
energy. This is demonstrated by plotting in figure 2 the moment .4,(/) of the flux in water 
as a function of energy in the range from 0.134 Mev down to 0.020 Mev for source energies of 
0.255 Mev and 1.0 Mev. The two curves have been normalized in such a manner that they 
have the same ordinate at 0.134 Mev. It can be seen that they lie quite close together; the 
area under the corresponding curve for 0.66-Mev source energy, which must lie between them, 
can be estimated with an accuracy of 5 percent. (For the purpose of this estimate, only multi- 
ply-scattered radiation from the 0.255-Mev source should be considered because a 0.66-Mey 
source will produce only multiply-scattered radiation below 0.134 Mev.) 

Computations were carried out for the energies and angles of incidence indicated in table 1. 

Singly-scattered radiation was calculated directly according to the equations given in 
appendix A. Some moment calculations were also made and are discussed in section 4. 

Multiply-seattered radiation was treated by the moment-method according to the function- 
fitting procedure outlined in section 2. It is desirable that one can fit the flux for as many 
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TABLE | Energies and obliquity angles for which the energy flux was calculated 


Asterisks indicate the values for which the multiply-scattered flux was determined by the method of function-fitting with the use of the weight 
function ¢(2z 5/2(2 


energies and source obliquities as possible with the use of the same weight function. This not 
only inspires confidence that the choice of the weight function is correct, but more importantly, 
it makes possible smooth interpolation. A considerable number of weight functions were 
tried, and the results obtained with them are described in section 4. 

Preference was finally given to the function 


go z)= 2°!/* KK. o(2) = L502), (13) 


where As, is a Bessel function of the second kind with imaginary argument. This weight 
function was superior to other functions considered in that its region of applicability included 
most of the energies and obliquity angles of interest (see table 1). In the relatively rare 
instances of nonapplicability, interpolation with the aid of results obtained with some of the 
other weight functions was successfully used. In any case, it is shown in section 4.2 that the 
results obtained with different weight functions were generally in excellent agreement. This 
provides pragmatic justification for the moment-fitting technique.‘ 

The results of the sample calculations for a 0.66-Mev plane monodirectional oblique source 
in water are illustrated in figures 3, 4, and 5. Figure 3 (a, b, ¢, d) shows energy spectra of 
the flux of scattered radiation at various distances from the source plane for obliquity angles 
a=0°, 30°, 60°, and 90°. Figure 4 (a, b) shows the dependence of the total scattered flux 


FO ( FE, 2)dE on the source obliquity and the distance from the source plane. 


When the penetration of radiation as a function of the obliquity angle is known, one can, 
by superposition, determine the penetration of radiation from plane sources with arbitrary 
angular distribution. An example of such an application is the calculation of the ratio 


| cos «| dEF(E,z) 


r(cos a’, 2 : — . (14) 


| deose | “dEF(E,2) 


This is the fraction of the scattered integrated energy flux from an isotropic 0.66-Mev plane 
source that is due to source radiation emitted within a cone of opening e’ around the normal 
to the source plane. This ratio is exhibited in figure 5 as a function of cos a’ for various values 
of 


‘One suspects that there may be a deeper-lying reason for the wide applicability of Bessel functions of the type used in eq (13 According 

) Spencer (private communication) such functions have also proved themselves useful for fitting flux moments for a point monodirectional source. 

lhe flux is expressed as the flow of energy, in Mey per Mev of spectral energy per unit time, through a unit area of surface perpendicular 

lireection of motion of the radiation, i. e., through an isotropic detector with an efficiency independent of the direction of the incident radi- 

ition Che normalization corresponds to a source emitting one photon with an energy of 0.66 Mev per unit time and unit area, Distances from 
ource plane are expressed in mean free pat! mfp) of the primary radiation 
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4. Checks and Comparisons 


4.1. An Experimental Comparison 





Experimental evidence with which the theory could be compared is unfortunately quite 
scarce, in fact is nonexistent for the exact conditions of the present problem. The only relevant 
published work seems to be that of Kirn, Kennedy, and Wyckoff [5], who measured the pene- 
tration of 0.66-Mev gamma radiation from oblique plane sources through slabs of concrete. 

A comparison will be made by considering experimental and theoretical buildup factors, 
The energy buildup factor-a characteristic parameter widely used in shielding work—is 
defined as the ratio 
Total energy flux (scattered plus unscattered 


Unscattered energy flux 


The fractional energy-transmission through a slab is the product of the buildup factor times 
the quantity exp (—2/cos a 

The experiment of Kirn et al. and the calculation differ in a number of respects: (a) The 
calculation for an infinite medium allows for the possibility that radiation can penetrate beyond 
the detector and then be backscattered into it. In the transmission experiment, however, with 
the detector behind a concrete slab and with no further backing behind the detector, such 
backscattering is impossible. The buildup factor for a slab is therefore smaller than that for 
an infinite medium. In a Monte Carlo calculation [6] it has been found that the ratio of the 
two buildup factors becomes nearly constant for penetrations greater than 4 mfp. (b) It is 
only approximately true that concrete may be considered a water-equivalent medium (by 
expressing all distances in mfp of the primary radiation). Actually, there is enough difference 
at low energies in the respective absorption coefficients (per electron) so that the buildup factors 
tor concrete are appreciably lower than those for water. By referring to the extensive tabula- 
tions of Goldstein and Wilkins [2] it can be found that this reduction factor is a very slowly 

irving function of the penetration distance (except in the vicinity of the source 
The reason for the approximate constancy of the buildup reductions due to effects (a 
b) lies in the fact that both are reductions of the low “CNerey end of the valima-ray spec- 


this part of the spectrum —as pointed out In section 3——1s an approximate equilibrium 


: trum, 1. e., rather independent of the depth of penetration. 
The calculation assumes a monoenergetic source spectrum. In the experiment the 
: vas physically large, and the source radiation was unfiltered, so that it presumably 
aun appreciable amount ol low-energ\ scattered radiation. Additional seattered 


could have been introduced by the collimating arrangement. 
The calculation is for the energy flux, whereas the experimental detector response is 
to the enere’s flux for energies down to 200 kev only, but rises somewhat more 
: t iower energies 


Hh shows semilogarithmic plots of the experimental and theoretical buildup factors 


sus oblique penetration distance 2/cos @ (in mfp , for obliquity angles a=0° and a=60 
seen that in both cases the experimental curves lie below the theoretical curves, and 
cos a) 5 the corresponding curves are almost parallel to each other, in conformity 
xpectation of a constant buildup-factor reduction. For medium and deep penetration 


refore correctly predicts the relative attenuation of the energy flux at different 
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4.2. Comparison of Different Weight Functions and Moment-Fitting Procedures 


In the absence of much experimental evidence that would provide a check on the theory, 
one is forced to rely heavily on tests of internal consistency. Hence a good deal of numerical 
experimentation has been accomplished, employing various weight functions to be used in the 
method of function-fitting, including also polynomial expansions whenever possible. 

Listed below are the relevant properties of a few functions that were found to be useful as 
weight functions in fitting the flux distribution from monodirectional plane oblique sources. 
These functions were selected because they have the proper asymptotic behavior and are either 
available in tabulated form or easily evaluated directly. In the formulas below, 2 will be under- 


stood to mean 


Bi, esel Function fs 


y(z)=L,(z)=2°A,(2), (15) 


Mo, dz gael) ad | ( 4 th ) T ( 4 +n-+p ): (16 








ov 
Lim L,(z (5) Pm 8154 





If p q- ,. (q ie 4 


L, (=) 51 ITY! 1 J9) “a . (18 
2 i=0 (q ayia! : 
(ie nerali ed Erponential Inte gral: 
. alt 
Yo / } é j 19 
(For p=1, this is the ordinary exponential integral). 
- f 
Lim /,(2) (20 
$n)! 
m (2] 
2n p 
krror Funetion: 
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yo ert(, ) | au ( dt 99) 
\7e 
Li 1—erf(, 2) “ 
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Table 2 gives a comparison of moment-fitting results obtained at an energy of 0.184 Mey 
for obliquity angles a=0° and 60°. Results for the symmetric and antisymmetric components 
are shown separately. For comparative purposes the normalization has been chosen so that 
the results obtained with the weight function Z;, are the standard of comparison and have the 
value unity. The over-all conclusion can be reached from an inspection of table 2 that the 
various methods are in good agreement: nonnegligible discrepancies occur only very close to 
or very far from the source. It is to be noted that at a=0° there is a good agreement with the 
well-established method of polynomial expansion into (= and )-polvnomials. Even a Gauss- 
ian weight function vields essentially the same results as the others at distances not too far 
from the source. The results of table 2 are typical and have been confirmed by other similar 
calculations. 

It should be mentioned that function-fitting could also be carried out with the use of a 
larger number of moments than four by representing the flux by a functional series of more 
than two terms. If the use of the additional moments makes little difference, the fitting 
procedure is presumably accurate. Such a check would be equivalent to establishing the 
validity of a polynomial expansion by examining its convergence. In the sample calculation 
the number of available moments unfortunately was too small to do this. 
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TABLE 2. Comparison of the results of moment filling with different weight functions 


Weight function (normal incidence: a=0 Weight function (oblique incidence: a= 
é 
0] 
I I I G hear, I I I I I I erf(y 
expansion 

1.00 0. UR 0). a7t 0. G18 1.06 0 1.15 1.04 1.00 0. O76 0. 961 1.30 1.12 1.09 

1.00 yu? Ist) “4h 100 0.5 0. UST 100 oo 44 OG4 1.03 oO Yu2Z 0. UUs 

Loo | 2 1.00 Lo 0. 478 l U36 0. Y62 100 Lo] 1.02 0. 920 “40 Y3yu 

2 100 0). WUS Lol 1.07 UTS 2 1s WUT 1 oOo 1.00 1.02 911 YUH ug 
‘ 1.00 7 0.449 0. GOS “74 { 1.08 1.02 Loo 0. US4 0.97 903 1.02 YHS 
_ 1 oo it) 1.01 1.12 O77 s 0. 404 0. O26 1 oo vn) .972 780 115 ao! 
100 i] 0.849 0.212 1.{ 16 1. 43 1.19 100 680 532 1. 20 1. 58 1™) 

Is4 f Is4 
! I Fi be I I I I I ! ! 
expans n 

1. 1.¢ 1O4 1. 51 0 1 ) 1.11 1.00 0. 920 0. 833 1. 46 1. 25 

9 1.00 1.1 109 1. 25 l 1.02 1.02 Loo 7s M51 0. 92 1.06 
1 oO 1 ou 110 1.04 2 1.12 0). Ut 1.00 108 1.02 410 O. 877 

j Loo 1 ol 0.433 1.11 1 1. ¢ 1.02 1.00 0.978 0. 950 1. 06 1.19 
* Loo 0. 82 1.11 O.S81I s 0. 447 0, 975 1.00 1.06 107 0. SSI 1.32 
lt 1.00 iol 0, 695 S46 If 1. 47 1.20 1. 00 0, 823 0. 698 1.63 1.78 

Is4 O.184 


4.3. Other Comparisons 


The penetration of radiation from a plane isotropic source can be determined through the 
superposition of results for plane monodirectional oblique SsOUrCceSs Good agreement has been 
found between the results of such a superposition and a direct calculation based on a polynomial 
expansion Wil h the use of the appropriate plane isotropic moments The ratio of the respective 
results for the total scattered energy flux at various distances from the source is shown in table 3. 

Asa further cheek, moment calculations were also carried out for singly-scattered radiation, 
which could be compared with the results of direct evaluations of the formulas of appendix A. 
Table 3 also contains sample results of such calculations; again, the agreement is good, It 
should be mentioned that the direct evaluation of the flux of singly -scattered radiation was 
rather laborious because of the necessary numerical integration. It would have been more 
expeditious, and would have involved only a minor loss of accuracy, to use moment-fitting 
making use of weight functions suggested by the exact analvtical expressions or physical 


considerations 


TABLE 3 Sample calculations to test the accuracy of the moment-fitling procedure 





1 F ( | 
| 
| 

i ih SP #. | 
; 74 } | 
0. 471 1.14 | | 
s “1s ] | 
res sf 8 iftered ff oma monodirectional plane ofl € SOUTCE 
A: Moment iation/direct éalculation (@=0°, E=0.325 Mev; weight function J 
B: ilatior t culatior r= E=0.215 M \V t function I 
~ ad | ty 
K) at 
( - t jue source result lirect calculation base 1 on plane otropic moments, 
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5. Appendix A. Singly-Scattered Radiation 


Consider a photon emitted by a plane source at z=0, with energy / and obliquity a. 


\ssume that the photon travels to a depth s without interaction, undergoes a Compton scatter- 
tween s and s+ds that changes its energy from £/) to a value in the range (2, E+dE), 
sses the c-plane with obliquity @ without further interaction. The probability for 
S venev is 
PAK Us K(k, E = ( whys uF): dy [(z—s) cos wldEds, (Al) 
COS a@ ( COS a cos@w 4) 


is the narrow-beam linear total-absorption coefficient, A(/, /) is the Klein-Nishina 
tial coetlicient for Compton scattering with an energy change from &, to /, and U’(: 


: init step funetion ((’°=1, if ->0; U=0, if -<0). According to the law of Compton 
mec mc ' Mic mc 
COS w=COS a COS{ l- » +a, SIN a sin ( l—— + ) cos D, (A2) 
( hk Ko ) eB Eo . 


here me? is the electron rest energy, and ¢ is an azimuthal angle distributed uniformly be- 
O and 27 

In order to obtain the energy flux /’, (/,2) of singly-scattered radiation, one must multiply 

eq (Al) by &, integrate over all values of s between 0 and © (assuming that O0<a<7/2), 

multiply the result by see w *, and average over all possible azimuthal angles ¢. With distances 

expressed in mean free paths of the primary radiation, the final result for the singly-seattered 


tlux can be written 


FO (BE, 2) Fo) caseon ay] |“ dell(z, Ea), AB 
Ll 1D COS @ ZT. 
where 
H(z, FE. a) } - : c 1) (cos w)—U COS @ > >U 
/ . . { - ) { 
ae =, COS w >< )). A4) 
] vi I | 
t ee ; 
COS @ 7 Ik COS W 


The integral in eq (A3) has been evaluated numerically, except W hen a=0 ; in this case 


cos w=cos(1— C4), AS 


and the integration over ¢ is trivial. 


6. Appendix B. Calculation of the Spatial Flux Moments 


The following brief discussion of the moment equations is presented to make this paper 
self-contained. It will also bring out some minor points in which the procedure followed here 
differs in details from that given in the basic reference (1). These are the use of an energy flux 
instead of a number flux, the separation of unseattered and singly-scattered flux, and the 
extension to plane sources with arbitrary angular distributions with cylindrical symmetry. 


? Without the factor sec @ one would obtain an energ urrent across the 2-plane instead of the energy flux through surface perpendicular to 
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Let F(z, cos 6, F) be the energy flux of photons that have been scattered & times. In terms 
of the F;’s, the transport equation has the form 





F , 
h cos A uk, SS’ > KE. Cos A) 
O 
OF, : *E (nee 
= cos @- uF | dE a5 K(k EF). 
; : ° > (B1) 
“2m “1 l me- . me* 
dy’ d cos 0’ 6 ( |——.. +; — cos 8 cos 0’ — 
|, , | l Qn ke iD 
sin @ sin 6’ cos e’ ) F, ' (k >0) 
4 
where 6 is the Dirae delta function, and S is the source function. 
We shall consider a set of elementary source functions 
: 2h+1 (B2) 
S,=ho6(h—Eo)P, (cos 0)6(z £==6,1,3....., 


Lor 


where PP, Is a Legendre polynomial. These elementary sources have no direct physical 
meaning except when L=0 but are convenient mathematical fictions. 

We shall indicate the source type from which a given flux is derived by an index L. 

We expand F,, into a Legendre series in the angular variable and take spatial moments of 
the Legendre expansion coefficients; 1. e., we introduce the transform 


Gn LE) | =| d cos 02"P; (cos 0) Fy, (2, cos 0,F) (B3) 
and a similar transform for the source function 


Tn dz d cos 02"P, (cos @)S,(2, cos 0, 3 6 (FE — E29) 6,0 bi. (B4) 
T 


l 


It follows from the transport eq (B1) and sundry properties of Legendre polynomials, that the 


transforms G,;, obey the equation 


A+ p)Go.=T, (k—0) 
. = . , Ti hy > mec mec . ns 
A+ w)Guy i} UE yey KU BYP (a en) Gace (k>0), (BS) 


where A is an operator defined by the relation 


; ’ /+-] : / 
AGy, nlk&)=n ; TT aACL 1./+1,.F)+n 


2-4] , Spy Haak 1,/—1,E). (B6) 


The transform of the multiply scattered flux is found by summation 


Hn, E)= Gen l,E) (B7) 


It follows from eq (B5) and (B7) that 


: "a — eo ne~ . me - "s - . 
(A+ p)ll,(n LE) | dk po KE way 1 ( 1 + )H(nl,E + T#(n,l,E), (Bs) 
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where the source function 


T#(n JE i} dE! 5, KE" EYP, ( — + ) Gln LE"). (B9) 


The calculation thus proceeds as follows: One computes Gz and G,, with the use of eq (B5), 
then the source funetion according to eq (B9), and finally H;, according to eq (B7). The equa- 
tion for Go, is algebraic; all others are linear integral equations of the Volterra type, and are 
solved numerically on the NBS automatic computer (SEAC). The equations are interlinked, 
but the nature of the coupling operator A is such that they need not be solved simultaneously. 
Indicated below are the chains of equations that must be solved to calculate //,(n,/,2) o1 


G.,(nJ/.E); for the sake of brevity only the (n,/) values are shown. 


i 


00 02 O4 O06 
1] 13 1] 15 13 15 
22 20 24 22 20 24 22 24 
; | o | > ol 33 
12 1) 12 10) 12 10) 12 
a1 a1 a1 51 
ov 60 60 60 
L=0 L=2 L=4 L=6 
Ol O3 OS O7 
12 10 14 12 16 14 16 
2 21 20 23 21 29 23 25 
| 52 0) 4 2 4) 4 b2 34 
1 tl 15 11 1 tI 3) 
52 50 52 50 52 50 52 
61 61 61 61 
70 70 70 70 
I I I ) L=§5 L=7 


For example, in order to compute //,(4, 0, /), one must first calculate in turn the quantities 


H1,(0, 0, BF), Hh, 1, F), Hy(2, 2, FD, and (4(3, 1, F 


It ean be shown that 


Hn LL.E)= Gn 1. E)=0 B10 


unless 7 L—/ is a nonnegative even integer. 


Now consider a monoenergetic plane source with arbitrary distribution, expanded into a 
Legendre Series 


. ; ~ 2L4+1 
‘S » COS Ak iD 0 ip ip i) >> Af, COS ff). Bil 
i tor 
It follows from the linearity of the transport equation that the moments of the flux for such a 
source can be obtained through the addition of moments for “elementary” sources. Taking 


into account B10), one finds that 
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(E)= >) Ao, He,(2n,0,F 
L=0 


. B12 


(E)=D>3 Aon, Hoy 4)(2n+1,0,E). 


I 





4 


it is remarkable that in order to determine the first n even (odd) moments, one need specify 
the source only by its first n even (odd) expansion coeflicients Az. 
Finally, we note that for a monodirectional source with angular part 6 (cos @—cos @), 


Az P, (COS @). (B13) 


7. Appendix C. Details of the Function-Fitting Procedure 


For the solution of the moment eqs (9a, b) it is convenient to adopt the following 
notation: 


A,=a,B a C,;=7,D = } 
Mi Mt, 
((* ) )) 
R mo M, ; nis wl D m, M,. |" (C2 
I~) m, M,' m, M, ° J 
j=1,2) (j=1,2) 


Equations (9a, b) can then be written in the form 


on} en 1} 1 +%=1 ) 


a8, a2B2= |] ¥16; + Y262=1 
he: A (C‘4) 





a) Pj > A2P5—= 7:6;7-¥265= 8 
a 
01; 9) > Abe P5 ty Y)0 Y205 ty 


Equations (C3) and (C4) being formally the same, further discussion is confined to (C3). 
A solution can be obtained by a method given by Chandrasekhar [7]. Let 


f g~ ‘ SN t ‘ 
Ww C5 re) (6) 
XN ] ‘ | 


It follows from (C3) that 8; and 8 are the roots of the quadratic equation 


B°+- w, B+ wo=0, (C7) 
and that 


B, and B, are real-valued if and only if 8, and g) are real and positive. This requires 
that 


w, <0 
Wy 2 U . (C10) 
Ww bw 0 


Equations (C5), (C6), and (C10) imply the following conditions on s and ?: 
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IfOSsS1,0StS 1, the allowed region lies between the parabola s?=t, and the curve whose 


parametric representation is 


1 - a Oe 
\ 3+ 6u—u*—|(3+-6u—u*)°—64u]'/" 
S (C11) 
t US 
(O<u<1) 


If s>1 and ¢>1, the allowed region lies between the line s=1 and the parabola s?=t. 


The author expresses his gratitude to L. V. Spencer and U. Fano for helpful discussions and 
suggestions. He also is indebted to J. Doggett and J. Hubbell for their help with the numerical 


calculations 
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A Re-entrant Cavity for Measurement of Complex 
Permeability in the Very-High-Frequency Region 
R. D. Harrington, R. C. Powell, and P. H. Haas 


\ re-entrant cavity, constructed for measuring the initial complex permeability of 


magnetic toroids in the very-high-frequency 


region, is described. Many of the electrical 


and mechanical problems associated with the development of an instrument of this type are 


considered. 


The general theory of permeability measurements, 
is discussed, and the associated working equations are derived. 


using a re-entrant cavity, 
Calibration techniques 


and results of measurements on several samples are also given. 


1. Introduction 


The increasing use of magnetic cores for high-fre- 
quency inductors and tuning applications has created 
a demand for information concerning the magnetic 
properties of materials in the very-high-frequency 
region. Research and production techniques are thus 
in need of a rapid and reasonably accurate method 
for measuring the initial complex permeability of 
toroidal samples in the region of 100 Me. 

Corresponding measurement problems in the fre- 
quency region below 30 Me were solved at the 
National Bureau of Standards through the develop- 
ment of the radio-frequency permeameter |1].!. Meth- 
ods used in the analysis of the permeameter are not 
readily extended into the higher frequency ranges 
because it is well known that an accurate measure- 
ment of radio-frequency permeability in the VHF 
region requires distributed parameter systems. These 
techniques usually involve impedance measurements 
with coaxial lines or cavities containing a sample 
whose dimensions are small compared to the wave- 
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Figure 1. Schematic drawing of re-entrant cavity for deter- 
mining radio-frequency permeability of magnetic toroids. 


Figures in brackets indicate the literature references at the end of this paper. 


length within it, and which is placed in a region of 
high magnetic field. When consideration is given 
to the study of permeability samples of very low loss, 
a cavity measurement is most readily utilized because 
of the attainable high Q of the instrument. For this 
reason, an investigation was undertaken for the 
development of a re-entrant resonant cavity for the 
measurement of initial complex permeability of 
toroids in the frequency range of 60 to 180 Me. 
A schematic drawing of a re-entrant cavity of the 
type developed under this research program is shown 
in figure 1. 

Cavities of this type have been constructed by 
other investigators [2, 3, 4], and as such the basic 
principles of these measurements are fairly well 
known. The design and analysis of the instrument 
developed at the National Bureau of Standards 
differs somewhat from these previous models. A 
description of the cavity and a brief discussion of 
the measurement methods and calibration techniques 
used in the determination of permeability with this 
instrument are given. Many of the electrical and 
mechanical problems that arise in the construction 
of a re-entrant cavity for magnetic measurements are 
discussed. 


2. Description of Instrument 


A re-entrant cavity is, essentially, a shorted section 
of coaxial transmission line with a lumped capaci- 
tance at one end, as can be seen in figure 1. The 
capacitor for the Bureau instrument consists of two 
parallel capacitor plates whose capacitance or sep- 
aration distance can be adjusted by means of the 
vernier micrometer drive shown near the top of 
figures 1 and 2. The complex permeability of a 
sample is determined from a measurement of the 
variation in @ and the change in the resonant capac- 
itv of the cavity when a small toroidal sample is 
placed at the bottom of the instrument in a region 
of high magnetic field. It should be noted that 
measurements can also be obtained by holding the 
capacity fixed and determining the change in 
resonant frequency upon insertion of a sample [2]. 
However, this variation in frequency is quite small, 
and its measurement requires elaborate frequency- 
measuring equipment, such as counters or hetero- 
dyne systems. As an instrument was desired that 


129 








a fairly accurate permeability measure- 
ment for laboratory as well as production and 
col trol purposes, emphasis was placed on 
capacity-variation techniques that require only an 
oscillator and detector as associated equipment, 

If the cavity is to cover a reasonably large fre- 
queney range, the tuning capacitance must be 
capable of a variation of the order of 10 or 20 to 1. 
to increase the frequeney coverage further, 
the mstrument was constructed in two lengths. as 
shown in figure 2. The longer instrument is used in 
the 60- to 110-Me region, whereas the short cavity 
110- to 180-Me region. The flange 
removal 


sections of 


Would ive 


qualits 


However. 


resonates in the 
in the middle of the cavity facilitates the 
and interchange of the two bottom 
different lengths. This makes it necessary to con- 
struct only one top unit, consisting of the capacitor 
plates and vernier unit. 

In order to avoid making the unit excessively long, 
the maximum capacitance must be large if low 
frequencies are to be obtained. This capacitance 
wes obtained on the National Bureau of Standards 
instrument by using capacitor plates 1.5 in. in diam- 
eter. Placing 9 toroid at the bottom of the cavity 
would require slipping r it over these | lates. This is 
undesirable as few representative cores have inside 
dimensions of 1.5 in. To insert the core through 





Ficure 2 Re-entrant cavities for very-high-frequency permea- 
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the bottom of the « ‘avity would require removing 
some type of shorting plug in the bottom section of 
the instrument, which, in view of the need of a high 
reproducible (, is impractical. That region of the 
cavity is, of course, a current maximum, and ch unges 
of resistance at the bottom cause large changes in 
cavity losses. 

This problem was solved by using a hollow center 
conductor through which a copper rod could slide. 
A gap in the hollow cong was placed at a con- 
siderable distance above the bottom of the cavity 
and the center copper rod completed the circuit by 
making contact with beryllium- -copper fingers on the 
inside of the hollow conductor. This opening in the 
center conductor may be seen the lower window 
of the two cavities in figure 2 and, in more detail, 
in figure 3. A core may be inserted by removing 
the cap at the bottom of the cavity and withdrawing 
the copper rod. The sample is then placed on a 
small plastic platform held by nylon threads and 





Figure 3 Lows rindow of re-entrant cavity sho na tech- 
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lowe! d to the bottom of the cavity. This platform 
and « typical sample can be seen in figure 3. This 
method of sample insertion thus allows the center 
conductor to have a much smaller diameter than the 
condenser plates, and, of course, a more reproducible 
0 is obtained. 

The capacitor plates can be seen in the top window 
of the cavities in figure 2. The removable window 
readil\ facilitates the cleaning and Inspection of the 
condenser as the capacitance is quite sensitive to 
dust particles that mia collect on the surfaces of the 
plates The lower condenser plate is held in place 
by two dielectric slabs, and a bellows connection is 
made with the top plate in order to avoid any sliding 
contacts. The capacitor plates are lapped to within 
0.0001 in. total indicator reading. The micrometer 
wheel can be read to 0.0001 in. directly and a vernier 
to 0.00001 in. This accurate dimensioning is neces- 
sary as small capacitance differences are used in the 
calculations. The calibration of the capacitor will 
be explained later. 

Coupling to the cavity is obtained through two 
adjustable current loops. The input loop is located 
near the top of the cavity, and the pickup loop for 
the detector Is placed a short distance above the 
sample, at the bottom of the instrument. Coupling 
must be loose in order that the input voltage and 
frequeney remain unaffected as the cavity is tuned 
throughout the resonance curve. A lyo-amp fuse is 
used as a bolometer element in order to obtain a 
better square law characteristic than is possible with 
a crystal, 

A typical equipment set up is shown in figure 4 
A signal modulated with a 1,000-ceps signal is passed 
through a filter and three-stub tuner. The input 
voltage is monitored with a vacuum tube voltmeter 
in order to assure loose coupling. The output is fed 
directly to a bolometer and amplifier. The fre- 
queney counter shown in the figure is used for eali- 
bration of the 0 of the eay ity. Once the () ealibra- 
tion has been completed, the counter is no longer 
essential for a measurement and is removed from 
the circuit The calibration details are given later. 


3. Analysis 


It has already been stated that a permeability 
measurement is based on the change in the inductive 
properties of the cavity when a sample is placed in 
a region of high magnetic field. The method is 
essentially 2 perturbation technique that enables one 
to disregard the effeet of the dielectric constant of 
the specimen 2s long as the thickness of the sample 
is small compared to the wavelength within it. It 
should be noted that this situation is satisfied in the 
VHF region by many of the samples manufactured 
commercially, In all cases, however, one must still 
give consideration to the problems that may arise 
with certain samples because of dimensional res- 
onance [5]. 

Consider the inner and outer radi of the internal 
volume of the coaxial section of the re-entrant cavity 
to be A and B, respectively. Let a, 6, and h, respec- 
tively, be the corresponding inner radius, outer 


radius, and height of the sample being considered 
If the thickness of the sample is small compared to 
the wavelength within it, and if the sample is placed 
against the shorted end of the Cavity in a region ol 
high magnetic field, one ay consider the electric 
energy stored in the sample as negligible compared 
to the stored magnetic energy, and write for the 
latter quantity 


, ve ' 
[ Hi 5 \K 1)L, f,\* 


where L, is the geometrical or equivalent air induc 
tance of a sample of permeability uv’ and is given by 
L,=0.4606 h logy, (6/2) * 107° henries, where / is in 
meters. /, is the current flowing in the vicinity of 
the sample and is essentially constant over the 
distance h. 

If we similarly consider the annular region of the 
empty cavity of radi A and / and height A/ located 
at the shorted end of the instrument, where A/ is 
sufficiently small, so that the electric energy in this 
region is also negligibly small, then one may write 
for the stored magnetic energy for this case 


yy = 

“ 22r A 
as the inductance per unit length of a uniform coaxial 
transmission line is given by 


where uo is the permeability of free space. 

It can be seen that there will be some value of 
AM/ for which 1 y,=ly2. By equating eq (1) and (2), 
this value of A/ is found to be 

; rl, 
Al=(p"—1) ; (3) 
] B 
olit 
_ 


This equation essentially relates the effect of 
inserting a sample of permeability yu’ into a region 
of high magnetic field in the cavity to a virtual 
change in length, A/, of the instrument where both 
the sample thickness and A/ are small compared to 
the wavelength being considered. 

This change may be transformed into a corre- 
sponding change in the lumped capacitance, C, 
located a distance, /, from the short by means of 
the approximate relation describing a short-circuited 
high-(Y transmission line with a lumped capacitance 
at one end, which is given by 


. | 
Z,, tan Bl Mad (4) 


where 
Zy—characteristic impedance of the cavity con- 
sidered as a coaxial transmission line 
B=27r/wavelength 


w=angular frequency. 
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Holding the frequency constant and differentiat- 


ng, we obtain 
who . 
d=— sin? BUC. 
Substituting the value of A/ found in eq (3) for dl 
In the above equation, we gel 
Ink wh, ° ’ 
=— sin® BIC. 
B re 
uo In 


1 


Solving for u’—1 and using the expression for the 
characteristic impedance of a coaxial transmission 
line Z, | /2m)y uo/€o In (8/4), where € is the per- 
mittivity of free space, we obtain 


__ de ; 
v hep Ly’ 9) 


where 


k.=(Z? sin? Bl). 


The above relation fives the real permeability in 
terms of dC, which is the change in capacity necessary 
to reresonate the cay ity upon insertion of the sample. 
The evaluation of fk, is explained later. 

The working equation for the measurement of the 
imaginary part of the permeability (u’’) may be 
tained directly from an analysis of the energy asso- 
ciated with the magnetic field in the cavity. If the 
usual relations for the magnetic field energy stored, 
as well as that dissipated per evele in both the cay ity 
and the sample are substituted in the corresponding 
expression for the ( of the loaded cavity, a relation 
is obtained that readily leads to the expression 


See ; 0) -,/@Q 
mM | 1)--*+ ‘( —1 ) , (6) 
‘a Gr La\Q, 


Z [3 sin 2~l | gl 


) ‘ 
a W 


ob- 


where 


and 
(/-=quality factor of empty cavity 


(»=quality factor of cavity plus sample. 


For a high-(V loaded cavity, the ratio of (’s may 
be replaced by the corresponding voltage ratio 
Vr= V, Vz, which is the ratio of the voltage out- 
put of the empty cavity to the output of the loaded 
cavity. Thus, the equation for u’’ is given by 


‘7 | 


‘ = 
y,| 1) Vre L (V7 1} (S) 


We therefore see that once yu’ 1 has been deter- 
mined for a specimen, the only additional informa- 
tion necessary for calculating the magnetic loss is 
the voltage ratio Vr, the @ of the unloaded cavity, 


See appendix 





and the value of the quantity 4, which is discussed 
shortly. 


4. Calibration of Cavity 


The re-entrant cavity is calibrated directly ugainst 
the primary standard of Radio-Frequency Per. 
meability at the National Bureau of Standards. 
which is a precision variable-length coaxial line. 
There are essentially four quantities on the re-entrant 
cavity that require calibration. 

1. The vernier micrometer drive, which must be 
calibrated in terms of capacitance. 

2. The quantity k in eq (5), which determines the 
real part of the permeability. 

3. The value of k; in eq (8), which is the working 
equation for the magnetic loss measurement. 

4. The Q of the unloaded cavity. 

It is not necessary to have a calibration curve for 
@ if a frequency counter or similar instrument is 
readily available. However, as previously men- 
tioned, it is preferred that measurements be obtained 
without relying on frequency-measuring — instru- 
ments. 

Consider first the capacitance calibration. The 
vernier micrometer could be evaluated in terms of 
the calculated capacitance between the condenser 
plates. However, the effects of fringing capacitance 
must not be neglected. The capacitor is thus cali- 
brated at 1,000 eps, using extremely sensitive 
capacitance standards. This is done by removing 
the lower door in the cavity and opening the center 
conductor circuit by withdrawing the removable 
copper rod. The lower capacitor plate is now iso- 
lated from the body of the cavity and held in place 
by dielectric disks. The lower end of that portion 
of the center conductor in contact with the bottom 
capacitor plate may be seen in the lower window of 
the cavity in figures 2 and 3. An external wire 
connection may be made with the lower capacitor 
plate at this point. The upper capacitor plate is in 
contact with the body of the cavity. It is thus 
possible to connect the variable capacitor on the 
cavity in parallel with a precision variable capacitance 
standard. This parallel arrangement is _ placed 
across one arm of a capacitance bridge, and incre- 
ments of capacitance in the cavity capacitor are 
compared directly with corresponding increments in 
the standard capacitor. The calculated capacitance 
of the variable capacitor in the cavity for a plate 
separation larger than that usually encountered in 
measurements is taken as a reference level. The 
exact location of this reference point is somewhat 
arbitrary as one is concerned only with differences 
in capacitance in this instrument. The above data 
enables a curve of capacitance versus vernier dis- 
tance to be obtained; however, the capacitance 
covers a rather large range, which makes it difficult 
to obtain a precise graph. This difficulty is avoided 
by plotting only corrections to the calculated capaci- 
tance, which is obtained from the formula 


132 





wh 
d } 


tal 
’ 


alr 


ca 
tel 
(Z 
an 
ca 
al 


to 
Ee 
Lo 
of 
Ve 
al 
of 
ol 
pr 
0), 
Sa 
di 
di 


or 


gi 


pe 





issed 


1iNst 
Per. 
irds. 
line. 
rant 


t be 
: the 


king 


for 
t is 
len- 
ned 
tru- 


lhe 
; of 
ser 
hee 
ali- 
ive 
ing 
iter 
ble 
SO- 
uce 
ion 
om 
of 
ire 
lor 
in 
ius 
he 
ice 
ed 
re- 


ire 





where €9 is the permittivity of free space, and A and 
qd are the cross-sectional area and separation dis- 
tance, respectively, of the capacitor plates. 

The expressions for the quantities /, and ky have 
already been given, but are repeated below: 


k.=(Z," sin? Bl)! (9) 
, Bas 
Zo| sin 28 . at | 

ky = . (10) 


9,. 
me 


In order to evaluate these quantities, Zp» is first 
calculated from the usual expression for the charac- 
teristic impedance of a coaxial transmission line 
(Zo>=138 logy(B/A), where B and A are the outer 
and inner radii of the cavity). The value of the 
capacitance between the electrodes is then obtained 
at a low frequency, where the electrode separation 
is small and the resulting capacitance large enough 
to account for most of the cavity capacitance. 
Equation (4) is then used with the above values of 
Z, and C' te calculate a value of / the effective length 
of the cavity. A theoretical curve of k, and hk, 
versus frequency may now be obtained. The oper- 
ating eq (5) and (8) are solved for k, and /, in terms 
of vw’, w’’, dC, Q., and Vr. The wp’ and yp” are 
obtained from a measurement of a sample in the 
primary standard. The remaining quantities dC, 
(., and Vr are obtained from a measurement of the 
same sample, using the re-entrant cavity. This 
data is obtained at one frequency for each of the two 
different length cavities. The resulting experimental 
values of k; and /, are then used to adjust the the- 
oretical values. It can be seen that this procedure 
gives an effective value for the characteristic im- 
pedance of the cavity. 

The @ of a resonant circuit may be obtained from 
either of the relations 


fo 2C, 
i) 9 (11) 
Aj A€ 
where fy and Cy are the resonant frequency and 


total effective capacity, and Af and AC’ are corre- 
sponding 3-db band widths. In the case of a distrib- 
uted parameter system such as a cavity, the total 
effective capacitance cannot be calculated. Thus, 
the frequency relation is more meaningful. 

However, the i) of a cavity is quite sensitive to 
temperature variations, and a more accurate answer 
is obtained if the @ could be obtained at the time 
of each measurement instead of using one original 
calibration of Q versus frequency. In order that 
only capacitance variations be necessary for a 
measurement, it was decided to make an initial 
measurement of Af, fp, and AC. Equation (11) can 
then be used to calculate 2Cy, which may be plotted 
as a calibration curve. Any future determination 
of @ may then be obtained from the calibration curve 
and a measurement of AC. 
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two different samples. 
5. Measurements and Results 


The cavity was normalized against the primary 
standard by using a sample having a permeability 
in the order of 8, which is a typical value for several 
samples in this frequency range. A comparison of 
results obtained at 130 Me, using three samples 
measured in the primary standard and re-entrant 
cavity, are shown in table 1. 


Comparison of results, using a re-entrant cavity and 
primary standard at 130 Mc 


TABLE l. 


Primary standard Re-entrant cavity 


> imple 


, m1 ’ 
u tan 6 7 mv tan 6 
aa 
l 8. 102 0. 067 8.119 0, 065 
9 10. 645 oan 10. 768 . ORS 
5. O8 027 5.16 034 


A typical curve of permeability versus frequency, 
using the re-entrant cavity, is shown in figure 5 for 
a powdered-iron sample and a ferrite. 

The above data and curves show the instrument 
to be capable of good precision and fairly accurate 
results. This conclusion, along with the simple 
measurement procedure, gives strong support to the 
use of a re-entrant cavity for very-high-frequency 
permeability measurements. 
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7. Appendix. Derivation of Loss Equation | Substituting the above values in eq (12), we obtain 


Defi ot terms l B ~ I 
; (Ing )A In 


ime of cavity including space occupied by i a 1 
aTripoie (2, 7 , 
f l | = ~ bh 
me of sample In K- fu In—h 
Feta tiie te whiliais Wann Q.8 A Qa 
fal energy stored I space occupied t LTripoie 
dissipated per evcle in cavitv without Solving for »” a’ /Q,, we gel 
SaAnipie 
] energy dissipated per cycle in space occupied by | B . 
sample of volume 1 n 7 n 
LB nner radius, outer radius, and length of cavity ” AA , AK 
, , . Lu Tih l -_ \ 
r radius, outer radius, and height of sarmmple , Bh 4 ), i), () Bh h 
ength measured along axis of cavity or sample In ? ' ? In 
radial distance measured in cavity or sample a a 
phase constant in cavity 
0. = of empty cavity + , 
‘) O of sampk In B 
), 0 of cavity plus sample, ” l , (J. ' I A < td. 
_ u y’—1) **4 7 K(f —-1)) (a7 
Che ofl the loaded eay itv may be written 4 (, Bh In oe 
ad 
ke ks ro ~ 
VO,=2F = —° 12 
Ke. T Es, Now tl 
Let B= wy Moe Mm 
I] te 
Il cos bf, _ . : we 
/ The geometrical inductance of the sample Is 4 
re 


where //7 is the magnetic field at a distance ¢ from the wo [,,. 4 
T=} in? |) - 
2r a 


shorted end of the cavity. 


If we assume the thickness of a sample is small - 
compared to the wavelength within it, then we may | The characteristic impedance of the cavity is _ 
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Fundamental Factors Controlling Electrical Resistivity 
in Vitreous Ternary Lead Silicates 


Simon W. Strauss,' Dwight G. Moore, William N. Harrison, and Lloyd E. Richards 


The resistivities of selected vitreous ternary 


lead silicates containing ions of groups 


I, II, I11, IV, and also the ions chromium, manganese, iron, cobalt, and nickel were measured 


in the range 200° to 500° C 


and at an applied direct-current field of 525 volts per centimeter. 


In addition, the resistivities of vitreous silica and of quartz, with optic axis parallel and 
perpendicular to the applied field, were measured over the same temperature range. The 
resulting data are given as log resistivity-composition isotherms for the alkalies and alkaline 


earths and as log resistivity-temperature curves for the other ions. 


Correlations between 


resistivity and bonding energy and between frequency factor and heat of activation are 


pres nted 
of the activation energy principle. 


1. Introduction 


Studies of the high-temperature electrical resis- 
tivities of ceramic coating materials have become 
important from the standpoint of both theory and 
technology. The current investigation was <n out- 
growth of a project begun in 1951. At that time a 
request had been received from Army Ordnance to 
develop a high temperature, clectrically insulating 
ceramic coating suitable for the replacement of the 
more strategic natural sheet mica in special electronic 
components which must operate at elevated tem- 
peratures A literature survey revealed that condue- 
tion in vitreous silicates was believed to be caused 
mostly by the migration of alkali ions. Using this 
information as well as some empirical considerations, 
a number of alkali-free, high-resistivity ceramic 
coatings were developed. A paper discussing some 
of the results was published [1] * and the data showed 
that ceramic coatings merit serious consideration in 
the field of high-temperature electrical insulation. 

Although the ceramic coatings that were developed 
and tested proved suitable for a specific application, 
itwas thought that a thorough study of fundamental 
factors controlling electrical resistivity in’ simple 
systems should be undertaken, to provide 
knowledge that could be applied not only in the 
design of ceramic high-tem- 
perature insulation, but also in the elucidation of the 
conduction mechanism With these objectives in 
mind, a fundamental study was begun in 1952 under 
the sponsorship of the Office of Ordnance Research. 
This report covers the first phase of the study. 


lass 


coatings for use as 


2. Review of Conduction Theory in Glasses 


Investigations of the electrical properties of vitreous 
silicates were made as early as 1854 [2,3]. Seddon, 
lippet, and Turner [4] have summarized the status 


! ddress, Naval Research Laboratory, W ishineton, D.C. 
? Figures in brackets indicate the literature references at end of this paper 


The resistivities of the alkali lead silicates were interpreted from the standpoint 


in this field of study up to 1982. More recent ideas 
regarding the conduction mechanism have been pre- 
sented in articles by Morey [5], Murphy and Mor- 
gan [6], Stevels [7], Stanworth [8], and Eitel [9]. 

It is generally agreed by these investigators that 
conduction in glasses is ionic. It has been estab- 
lished that the electrical conductivity of a soda-lime- 
silica glass, at room temperature and up to the strain 
point,’ is essentially the result of transport of sodium 
which are believed to migrate through the 
interstitial positions in the silicate network. In this 
temperature region the mobility of the calcium ion 
is considered to be negligible. The thermal energy 
is thought to be sufficient to detach a sodium ion 
fromitsinterstitial position and cause it to move to a 
neighboring position [10]. As long as no external 
field is applied, the probability of such migration is 
equal in all directions. However, in the presence of 
an imposed field the sodium ions will show a tendency 
to move with the direction of the field, thus causing 
an electric current. The sodium ion has either a 
lower charge or a smaller radius than the other 
network modifiers ‘ in a glass, with the exception of 
lithium; for this reason its movement through the 
lattice Ta be less restricted. Nevertheless, at suffi- 
ciently high temperatures the other network modi- 
fiers are said to contribute to the conductivity [7]. 
Relatively few systematic studies have heretofore 
been made on the conduction of simple glasses at 
temperatures exceeding 200° C. Moore and De- 
Silva [11] studied the electrical properties of ternary 
and quaternary alkali lime glasses at temperatures 
up to 320°C, and Gehlhoff and Thomas [12] measured 
electrical conduction of ternary silicates at 
temperatures up to 400° C. The purpose of the 
present investigation was to provide information 
which might contribute to a better understanding of 
the conduction process in silicate structures, 


1OnS, 


soda 


3 Lillie, J. Am. Ceram. Soc. [3] 37, 117 (1954), recommended the strain point to 
be defined as that temperature at which the data from the annealing point test, 
when plotted logarithmically against temperature, extrapolate to a viscosity of 
10'4.20 poises, 

4 Network modifier ions usually occupy the interstitial positions in the glassy 
structure. 
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3. Selection of Compositions 


One of the objectives of the current investigation 
was to study the composition dependency of elec- 
trical resistivity of ceramic coating materials. It 
appeared that a logical approach was to select a 
simple silicate system for initial experimentation, and 
then to increase the complexity of the system by the 
cradually inereasing partial replacement of one 
component ion at a time by other ions. In selecting 
a base COMPOsItlon, the svstem chosen should have 
a high solubility for various oxides and should vield 
workable Careful consideration of these 
and other problems involved, resulted in the selection 
of PbO:SiO, as the base composition. In a given 
composi. lon series, increasing amounts of lead oxide 


classes. 


were replaced by one of the other oxides. 

The structure of a glass may be modified 
siderably if the composition is changed over a wide 
range even if the same constituents are used 
throughout. To avoid such structural changes, the 
specimens were composed according to the formula 
rMeO:(1 r)PbO:;Si0,. The value of « ranged from 
0 to 0.4 for the alkali and alkaline earth ions and for 
iron, cobalt, and nickel and was fixed at 0.4 for the 
other ions. The 0.4 mole replacement was selected 
because this was the maximum concentration for 
class formation for some of the oxides, and volatiliza- 
tion could be kept to reasonable values. 
Equivalent replacement was used so as to keep a 
constant ratio of Si/O and thus avoid any structural 
changes which might be caused by changes in the 
Si/O ratio. 

Resistivity measurements were made on quartz 
and silica glass and also on complex lead silicates in 
which part of the lead oxide was replaced by oxides 
of lithium, sodium, potassium, rubidium, cesium, 
beryllium, magnesium, calcium, strontium, barium, 
zinc, cadmium, boron, aluminum, lanthanum, tita- 
nium, germanium, zirconium, cerium, thorium, chro- 
mium, manganese, iron, cobalt, and nickel. The 
compositions of these glasses are given in table 1. 


cCon- 


le SSCS 


4. Materials and Preparation of Specimens 


All glasses were prepared from reagent grade 
chemicals (oxides, nitrates, carbonates, or  ox- 
alates); the source of silica was quartz crystals 
pulverized to pass a US Standard No. 40 
Each raw batch was mixed intimately and ground 
in a mortar until all of the powdered material 
passed through a No. 40 sieve. The batches were 
melted in platinum crucibles placed in an electric 
furnace. Each melt was periodically stirred with a 
platinum rod to aid in the elimination of occluded 
gases and to produce a homogeneous glass. Test 
specimens were cast in a steel mold in the shape of a 
disk, 1's in. in diameter by about 4 in. thick. The 
disk, after casting, was transferred to a preheated 
ceramic mold and was then placed in an annealing 
furnace where it was allowed to cool overnight. 
After cooling, the disk was ground flat on a grinding 
lap using SiC No. 150 and No. 320 as abrasives, and 
the thicknesses measured to the nearest 0.0001 in. 


sieve. 





TABLE l. Ratch « om positions of ternary lead silicates 


nlain. 


ing elements of groups IT, IT, IIT, and IV, and chromium 


} 


manganese, iron, cobalt, and nickel 


Composition Amount of 


No Oxide added poten PhO 
Moles Moles 
l Lio 0.01 0.49 
2 Li2O 10 Ow 
; LieO a) i) 
‘ I ) 30 70 
5 LiLo Ww) mw 
NaoO 01 ug 
7 Nao Th un) 
s NaoO™ *") si) 
; Nao +) 70 
10 NaoO Ww) mw 
11 KO Ol ag 
12 Ko it “ 
13 Ko 1) si 
14 Ko +) 70 
15 Koo 1) “ 
It) Rhee? HW) mw 
7 CoO Ww) mw 
Is Be) Hi) “wo 10 
14 MeO Ol au 
2 MgO Ww ") 
21 Meo " s() 
yo MgO ( 
y 4 MeO w) ‘ 
24 Cad 1 "7 
ys Cad Ww u) 
y Cad i) si) 
rai Cad ‘4 | 
Js ( 4) Hi) ow) 
pt | srt) 1 "uy 
4) sr) lt ”) 
l sr »”) st) 
5 <r) ‘ 
srt) Hi) 4 
‘ Bal) ] "0 ( 
Bad l " 
Bal) a si) 
BaO ‘ 
ts Bad HW) mm) 
su Ano HW) “wm 
1) cdo Ww) “wo 1.‘ 
11 BO 133 m) ) 
42 ALO 133 mw 10 
i] l 4) 133 4 0 
4 rio a) i) 10 
45 CeO, wa am 1.0 
1 ZrO a mw 
ri Ce a | A) ) 
s Pho ‘a til ) 
ae ( ‘) 13. “ 
4 Fe) T 4 ) 
| Coo ") mw) 
y2 NiO HW) i) 
Vino 1) “wo 


Specimens were prepared for test by spraying a 
thin laver of “silver paint’’’ onto face. A 
template was used during the spraying to leave an 
exposed circular area of silver *% in. in diameter 
evenly centered on each face of the specimen. The 
paint was cured by air drying at room temperature 
for several hours. The resulting silver contacts 
were smooth and continuous before and after test 
Results of a previous investigation [1] had shown 
that resistivity of alkali borosilicate and alkali-free 
borosilicate ceramic coatings appeared to be inde- 
pendent of thickness. It was therefore assumed 
that the contact resistance was small in comparison 
to the resistance of the specimen. 


ve purpose n the ceram 


* A term ipplied to a preparation used for decorat 
ifter curing, it consists @ 


industry and more recently in making printed circuit 
nearly pure silver 
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Several of the glasses were analyzed chemically. 
The compositions were chosen to represent high, 
intermediate, and low resistivity glasses. A com- 
parison between the compositions calculated from the 


5. Chemical Analyses 


hatch materials used and the compositions as 
determined analytically is given in table 2. 
| 
TABLI 2 ( he mical ana s of fernary lead s rales prepa ed 
hy ti equivale nt substitution of beryl ium, manganese, cobalt 
and ne oxide for lead oride in PhO:;SiO 
We t Weight percer 
{ 
( Cal 
0) le i Ovxice we Fou 
BeO ) 13 PhO ‘ ( 
ZnO t 14.4 PhO 9 1 ; 
Nino 2s 12.1 PhO 60 ‘ 
CoO { 11.8 PbO " } 
Che ] analy 1M e ana ic laboratory he M ] 
Produc ) on under the per n of L. Bean. I he lica content W 


6. Test Equipment and Procedure 


A photograph of the equipment used for determin- 
ing the resistivity of the specimens is shown in 
figure 1. The furnace, A, had a heating element 
consisting of Nichrome wire wound on a ceramic 
Within the core were located two circular 
steel in diameter. An 
iron-constantan thermocouple was embedded in the 
top electrode, and the control thermocouple was 
embedded in the insulation. This second thermo- 
couple activated the temperature control unit, B. 
The temperature of the top electrode was determined 
by the use of a potentiometer, F. A recording 
potentiometer, E, was used with a preamplifier, D, 
to graph the current passing through the specimen, 
1O~” to 10 amp. An Ayrton shunt, C 
was used in the input to the preamplifier to extend 
the range of the recorder to 10 amp. The resist- 
ance of the preamplifier input was 1,000 ohms. 

The specimen to be tested was placed in the fur- 
that the stainless steel electrodes were in 
contact with the silvered area on either side of the 
coated specimen (figs. 1 and 2). The furnace was 
heated at the rate of about 10 deg C per min and the 
current through the specimen was measured at 10 
degC intervals. The voltage was applied during the 
entire heating period, polarity being reversed at least 
once each 20 deg C. All determinations were made on 
at least triplicate specimens, one at a time. The 


core. 


stainless electrodes s In. 


In ranges 


nace so 


current through the specimen was measured in 
amperes and the resistivity (ohm-em) calculated 
from the current, voltage, and dimensions of the 


specimen. When the voltage drop across the shunt 
exceeded 1 percent of the input voltage or 2 v, this 
voltage drop was taken into consideration in eal- 
culating the resistivity. Kdge effects were not 
included in caleulatine the effective area of the 
specimen because of the general inaccuracy in this 
type of measurement. The silvered area was 


as- 
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Kiet RE | T sting arrange ment for the measurement of 
electrical resistivity of ceramic coating materials at elevated 
te m perature Ss. 
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electrical resistivity of ceramic coating materials at elevated 


temperatures. 


sumed to be the effective area. In this type of 
measurement where the resistivity values are high, 
a 100-percent difference between duplicate deter- 
minations may be obtained. The precision was 
found to be, in general, of this order of magnitude 
or better. Above about 200° C the same values 
for the current were obtained with and without use 
of a guard ring. The surface conduction was there- 
fore assumed to be negligible above 200° compared 
to that of the volume. 

The effect of polarization was directly observable 
in the values of current before and after reversal of 
the direction of flow through the specimen. In a 
specimen containing 0.4 mole of sodium, where 
polarization effects must be expected to be high, a 
maximum difference of about 6 percent was observed 
on reversal at the highest temperatures above about 
300°C. This level of error introduces no observable 
error into the activation energies calculated from 
these data, and has therefore been ignored. 








7. Results and Discussion 
7.1. Silica Glass and Quartz 


A comparison between the resistivity of silica 
glass and quartz is given in figure 3. These experi- 
ments were performed, (a) to measure the resistivity 
of silica glass because it is one end member of the 
PbO :SiO, system, and (b) to compare the resistivity 
of two allotropic forms of silica in order to obtain 
information regarding the conduction mechanism. 

The observed conduction through both quartz and 
silica glass might well be due, in part at least, to the 
presence of traces of impurities since a trace of 
impurity would vield many potential conduction 
ions. For example, an impurity of one part sodium 
per billion of silica would be equivalent to about 10! 
sodium ions per gram of silica. The results of 
spectrochemical analyses (table 3) indicate that each 
gram of silica glass or quartz actually contained 10! 
to 10'° ions of various impurities. While these ions 
would seek the lowest available potential energy 
holes in the silicate structure, it is quite conceivable 
that at least some of them might still contribute 
towards electrical conduction, especially at higher 
temperatures. Warburg and Tegetmeier [13] re- 
ported that when a quartz specimen was electrolyzed, 
sodium amalgam being used as the anode, sodium 
migrated through the specimen in) amounts in 
accordance with Faraday’s law, while the welght of 
the specimen remained unchanged. These investi- 
gators concluded from this result and from addi- 
tional experiments that the conduction was due to 
the presence of sodium silicate in the quartz; their 
analysis showed 1 part in 2,300, hence the quartz 
was regarded as a very dilute solution of sodium 
silicate 

As shown in figure 3, quartz 
measured perpendicular to the optic aXis Is greater 
than that of the same material measured parallel to 
the optic axis, over the range 200° to 500° C Quartz 
has two crystalline modifications, @ and 8. For 
alpha quartz, the existence of tunnels parallel to the 
optic axis has been suggested in the literature, while 
similar tunnels are postulated to be absent in the 
directions perpendicular to the optic axis [14]. This 
tunnellike structure is supposedly responsible for the 
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ease of ion migration in quartz plates with optic axis 
parallel to the applied fiell [15]. This hypothesis is in 
agreement with the results plotted in figure 3 

From figure 3 it is also seen that the resistivity 
curve for quartz measured parallel to the optic axis 
does not follow the Rasch-Hinrichsen law, log 
R A/T B, over the entire temperature rang 
that was stucied This effeet was observed on fou 
independent determinations using four different 
specimens. The change inh slope of the temperature- 
resistivity curve indicates that there is a change in 
the mechanism of conduction in the quartz whose 
Optic axis 1s paralle| to the applied field when the 
temperature exceeds about 360° C 

The transition temperature of @ to 8 quartz is 
viven as 573°C. However, X-ray diffraction studies 
performed on quartz [16] indicate that a partial 
reversible transition from the @ to the 8 form may 
occur between room temperature and 573°C. If 
this apparent behavior of quartz ts further sub- 
stantiated by additional X-ray studies, the results 
may be helpful in explaining the observed chang: 
in slope in the temperature-resistivity curve 

It is not within the scope of this paper to discuss 
the conduction mechanism in silica in any greater 
detail. The results do, however, indicate that a mor 
critical study of the electrical properties of allo- 
tropic modifications of silica might vield tnportant 
information regarding the mechanism by which ions 
migrate through a structure composed of SiO, group- 
In addition such a study might throw further 
structural arrangements of the 


ings. 
light on the exact 
various modifications of silica. 


7.2. Alkali-Lead Silicates 


Figure 4 shows the resistivity-composition  180- 
therms at 200°, 300°, and 400° C that were ob- 
tained when lead oxide was replaced with equivalent 
amounts of lithium oxide, sodium oxide, or potas 
sium oxide. (‘rossover of eurves which are similar 
to the ones shown in this figure are reported in the 
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literature for resistivitv-composition isotherms as 
well as for viscosity-composition curves for alkali 
silicates [17]. It will be noted in figure 4, that re- 
placement up to about 0.3 mole at 200° C, 0.2 mole 
at 300° Cy, and 0.1 mole at 400° C vield glasses 
whose resistivities are somewhat vreater than the 
resistivity of the base composition. Such behavior 
suggests that additions of alkali at the expense of 
lead up to a certain optimum amount actually 
tightens the structure, thereby making migration of 
ions more difficult. A tightening of this type may 
be expected if the lead is acting as a network former, 
a role that it is reported to ussume at high concen- 
trations [IS]. As more alkali ions are substituted 
for lead, however, av pont should be reached where 
the presence of extra alkali ions would more than 
compensate for the tightening effect resulting from 
the removal of lead from the network positions. 
The resistivity should then decrease upon further 
addition of alkali ions and figure 4 shows that such 


i decrease does Occur 


7.3. Alkaline Earth and Transitional Element Lead 
Silicates 


Figure 5 shows the resistivity-composition iso- 
therms for alkaline earth lead silicates, while the re- 
sistivity-temperature curves for iron, cobalt, and 
nickel lead silicates are given in figure 6. 

In figure 5, there is a definite trend toward in- 
creased resistivity with increasing concentration of 
the alkaline earth oxide. A behavior of this type 
might indicate a general tightening of the strueture 
from the replacement of lead with the alkaline earth 
ion, making it more difficult for alkali ions, present 
as trace Impurities, to migrate through the structure. 
A second interpretation might be based on the pos- 
sibilitv that lead and the alkaline earth ions are 
themselves responsible for at least part of the con- 
duction. However, such an interpretation would be 


Figure 5 


eaiciim, 


Resistivil j-com position isotherms for magnesium, 
strontium, and harium lead silicates pre pared by 
MeO. Cat ), SsrO. or BaO for 


equivalent substitution of 


PbO in PbO:SiO 
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y ranges from 


meaningless until it can be demonstrated by some 
suitable experimental technique that ions other than 
alkalies are free to move in a rigid silicate network. 

The slopes on the curves in figure 6 show that 
while iron and cobalt tend to decrease the heat of 
activation when 0.3;and 0.4 mole of these oxides 
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are added to the base composition, nickel tends to 
increase the heat of activation at all concentrations 
studied. The addition of small amounts of iron and 
cobalt appears to have the same effect upon heat of 
activation as a small amount of nickel. One pos- 
sible reason for this behavior may be that at small 
concentrations these ions (also the nickel ion), which 
are about as small as a lithium ion, enter deep po- 
tential energy holes and cause a tightening of the 
structure. The nickel glasses were also much more 
difficult to prepare than the iron or cobalt glasses, 
in that they showed a greater tendency to devitrify. 


8. Interpretation of Resistivity Data in Terms 
of Energy Considerations 


In view of the finite conductivity found in the 
quartz used as the raw material, it seems probable 
that the major part of the conduction in all the 
glasses studied may have been due to alkali ions, 
either added or present as impurities. This 
sumption will be implicit in the following discussions 
which will apply to ternary lead silicates in which 
0.4 mole PbO in a base composition PbO:;SiO, was 
replaced by oxides the ions of which are given in 
figure 7. 


as- 


8.1. Electrostatic Bonding Energy Considerations 


Various investigations [19] have used charge and 
ionic radius to obtain a semiquantitative classifica- 
tion of constituents. Such expressions as 
er, er, Zelr, or Zelr? (Z, e, and r representing 
valence, charge of the electron, and tonic radius) 
have been used to express the relative magnitude of 
cation-OXxy ven bond energy in vitreous. silicates. 
When other ions were introduced into the base com- 
position one of the results may have been a tightening 
or loosening of the silicate framework. The extent 
of this tightening (or, loosening) of the structure 
was probably dependent, in part at least, upon the 
size, charge, and the position occupied by the eation 
added to the base composition, The electrostatic 
bonding energy was therefore used as a first approxi- 
mation to estimate the relative tightening (or loosen- 
ing) effect produced by introducing other ions into 
The equation for electro- 
static bonding energy [20] is A 7? Ryo: where 
electrostatic bonding energy, 7; and Z, are the 
valences of the cation and the oxygen ion respec- 
tively, eis the charge of the electron, and /?). is the 
This equation states that 


glass 


the base composit ion. 


cation to oxygen distance. 
the magnitude of the bonding energies Increases W ith 
the valences and decreases with the atomic separa- 
tion Figure 7’ 1s a plot of the logarithm of the 
resistivity at 200° C against the logarithm of the 
quantity (—e-*F) which is proportional to the elec- 
trostatic bonding energy hg Possible oxidation states 


Since all glasses water it is« urse | bl 
ma Iso contribute towards conductiot 

? Similar curves were obtained at 300, 400, and 500° C 

§ Crystal radii were used in the calculation of electrostatic | ercie 
It is realized that these ionic radii values are relative. In a glass, wt yrrie 
time described a i upercooled l jul i, al olute value for ionic tn 
expected t lifferent from ecrvstal radii Valuc Chan 
number n ilso affect the re3 of ionic radil, 








lead 
silicate and for ternary lead silicates ore pared hy the equivalent 
replacement of O4 PbO in the hase com position PbO:SiO 
hy each of the ions listed on the curves. ‘ 


Ficgure 7. Resistivity-bonding energy isotherm tor 


itic bonding energy, e=charge 


ind oxygen ion, respec 


In the bonding energy expression, E =electrost 


of the electron, Z; and Z» are the valences of the cation 
tively, and R cation-oxygen distance 


are indicated for iron and manganese to illustrate 
the variable valency of transitional element ions. 

It is interesting to note that the ions added to the 
base composition, seem to fall into three general 
classes. The first class contains ions of groups I, 
Il, and the ions Al and La, the second group con- 
tains tetravalent ions, and the third, less defined 
group, contains the ions chromium, manganese, iron, 
and cobalt. With the exception of iron and cobalt 
only the alkali ions added to the base composition 
decreased the resistivity. This observation appears 
to be consistent with the assumption that the major 
part of the conduction may be due to the alkali 
ions. 

Figure 8 gives a plot of upper and lower limits of 
heats of activation against composition for lead sili- 
eate and for lithium, sodium, and potassium lead 
silicates. The heats of activation were calculated 
from the slopes of curves obtained by plotting the 
logarithm of the experimentally determined resistiv- 
itv values against the reciprocal of the absolute tem- 
perature. The data indicate that there is an initial 
rise in heat of activation followed by a drop at the 0.4 
mole composition. An initial rise in heat of activa- 
tion is in accord with the suggestion made earlier 
that initial replacement of lead by alkali might vield 
a tightening of the structure. 

A plot of heat of activation, AII*, against the log- 
arithm of the frequency factor, log A, is given in 
figure 9. The quantities log A and AH * appear 
as parameters in the equation log p=log A4 
(AH*/2.3RT) where p is the resistivity. Although 
the seatter is relatively large and log A is a complex 
quantity, the decrease in log A observed with the 
alkali ions contrasted to the constaney or slight in- 
crease in log A for the other ions is consistent with 
the alkali ions alone providing additional charge 
carriers. 


140 





Fr 


cu 
lit 
po 


be 
th 


OX 


en 
stl 
tu 
su 
di 
im 
an 


In 
fo 
th 
of 


ho 


trate 


is. 

» the 
eral 
s I, 
con- 
‘ined 
iron. 
rbalt 
ition 
ears 
sajor 
Ikali 


is of 
sili- 
lead 
ated 
the 
stiv- 
em- 
itial 
().4 
iva- 
rlier 


ield 


log- 
1 in 
pear 
A+ 
ugh 
plex 
the 
in- 
vith 


irge 











16 | 
aa 
34 
» 
Fs 32 » 
= » 
Q » 
e . 
7 
* 3 
2 » 
~ e 
. ; 
$*"o8 
if 8 
} 
> 26 
4 . 
. 
24} 
’ 
- - . 
22} 
I t 
| 
20 Li J 
20 3 4c 
ALKALI OXIDE , MOLE 
FIGURE 38. Upper and lower limits of heat of activation 


composition data for lead silicate and for lithium, sodium, 
and potassium lead silicates pre pared by equivale nt substitu- 
fion of Li 4). NaQ, or kK () for PbO in PbOrSiO». 


Figure 10 shows heat of activation-ionic radius 
curves. For the alkali ions, AH* from 
lithium to potassium and then begins to increase from 
potassium to rubidium. The minimum or at least 
change in slope at potassium appears to be real and 
beyond experimental error. It can be explained in 
the following way. 

Before an alkali ion can migrate, the free energy 
of activation must be exceeded. This involves over- 
coming the chemical binding energy (breaking the 
eation-oxygen bond) and, in addition, sufficient 
energy (network repulsion energy term) must be 
supplied to the alkali ion so that it can pass through 
the interstitial hole. From the standpoint of the 
chemical binding energy term, the smaller the cation 
the more energy it will take to break the ecation- 
oxygen bond. On the other hand, the smaller the 
eation the smaller will be the network repulsion 
energy term required to move it into the next inter- 
stitial hole. Assuming that the approximate struc- 
ture of glass does not differ appreciably from a model 
such as is depicted by Ponecelet [21], the average 
diameter of the opening between interstices is approx- 
imately equal to that of an oxygen ion. Anderson 
and Stuart [22] in discussing the structure of glass 
state “Although it is true that the volume of the 
interstice may change when an ion is added, there are 
four doorways connecting the interstice containing 
the ion to adjacent interstices and the average size 
of the doorway through which the ion must pass will 
not change significantly.””* An ion smaller than 


decreases 


« possible that some of the openings between interstitial holes may 
he circumference of an oxygen ion. Even if this would be the case 
would still be generally the same the one presented. 
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FiGcurRE 10. Heat of activation—ionic radius curves for ternary 


lead silicates prepared by the equivalent replacement of O.4 
PbO in the base composition PbO:SiO, by each of the ions 
isted on the curves 


oxvgen should pass through the opening to the inter- 
stitial hole with considerable ease whereas ions larger 
than oxygen have to expand the entrance to the inter- 
stitial hole before they can pass through. Hence, 
in the latter case the network repulsion energy term 
may become a large fraction of the total energy re- 
quirement and the ion would have to surmount a 
higher energy barrier before it could migrate into the 
next interstice. The ions rubidium and cesium are 
larger than oxvgen therefore the network repulsion 
energy term for each of these cations may be ap- 
preciably larger than that for any of the first three 
alkali ions. This increase in the network repulsion 
energy term apparently more than compensates for 
the decrease in the chemical binding energy term 
with a net result that the total energy of activation is 
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greater for rubidium or cesium than for potassium. 
This explanation should be accepted with some reser- 
vation until additional experimental data are 
available 

The preceding discussions may be summarized as 
The data indicate that the tons added to 
the base composition fall into four groups: (a) 
alkali ions tend to decrease log A and to vary the 
heat of activation, with an indication of decreasing 
the heat of activation with increasing radius, and 


' 
follows 


with a possible minimum or a change in slope at 
potassium: (b) iron and cobalt do not appear to 


affect log A but decrease the heat of activation: (¢ 
mat affects neither log A nor the heat of 
activation: and (d) the other ions increase log A, 
and increase the heat of activation, with an indica- 


‘ 


ranese 


tion that within a group of isovalent ions the heat of 





Relative to the base composition, substitution for 
0.4 mole Pb* of alkali ions, iron, and cobalt decreased 
the resistivity, manganese did not affect the resis. 
tivity, and the other ions increased it. 


The authors thank J. D. Hoffman of the Polymer 
Structure Section and A. Franklin of the Porcelain 
and Pottery Section of the Bureau for their helpful 
discussions on certain aspects of the problem 
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Some Factors Affecting the Precision of Polarographic 
Half-Wave Potential Measurements 


John K. Taylor and Stanley W. Smith 


The experimental conditions required for the precise determination of polarographic 
half-wave potentials have been investigated and are discussed. The values of half-wave 
potentials vary with the characteristics of the dropping-mercury electrode in the case of 
reduction of metal ions to metal soluble in mercury Consequently, it is nee*ssary to specify 
the experimental condition whenever values for half-wave potentials are reported. With 
the 


roper precautions it is possible to attain reproducibility within +0.2 millivolt ir 
pro} | | | 


measurements 


1. Introduction R, 


The polarographic half-wave potential is of in- ey 


terest because of its use in the qualitative identifica- 
tion of electrode reactions. Of even greater interest, i 
however, is its thermodynamic significance and its 


application to the problem of determining the com- a 
eo— 


plexity constants for complex ions [1,2]. In most 


eases the half-wave potential data presented in the (<, ) 


literature have been precise to less than plus or minus 












































several millivolts. While this is sufficient for quali- ©) 

tative purposes, it is not of the required precision for 

the other applications The purpose of this investi- C 11 

ration Was to determine the precision with which it Q— 

e . ; 7. 

s possible to measure half-wave potentials and also t — Ge 
to study the factors affecting the measurements. . \ REF ~~ J 





As the purpose of the work was to investigate 
methods yy which precise values Tay be obtained 
and not to establish a precise half-wave potential PiGure 1 Manual arrangement for the measurement of half- 
spectrum, relatively few reducible materials were wave potentials. 
examined 


2. Experimental Measurements 


All half-wave potential measurements were made 


with i manual arrangement, the cireult diagram of ——} 














which is shown in figure | The voltage was applied AT \ ] —— No 
to the cell by means of the slide wire (S) of a Sargent | | “———_~ 

model XILT or model NXT polarograph. The poten- . | ANODE 

tial of the dropping-mercury electrode (DME) with | | [Tt Ne tee 4 
respect to the saturated calomel reference electrode en cs Ba mt + 

REF) was measured with a precision potentiometer | —_—e iif 

Pot The current was measured with an Ayrton a a er ee ee 
shunt (SH) and a lamp and seale galvanometer | | 7 1 et aoe 

G Galvanometer (G,), the null indicator on the 

potentiometer, was a “spotlight” type of instrument. | 

A system of B., Ro, Ry, Ry, and R; was included for 18 

















the calibration of the galvanometer (G,). This 
calibration current could be calculated from = the 
voltage drop across R, and its known value of 


resistance 

















The cell arrangement shown in figure 2 consisted 
of a 100-ml beaker fitted with a rubber stopper that 
contained inlet tubes for nitrogen for the removal 
of oxygen from the solution and holes for the access 





— Figure 2. Polarographic cell used in half-wave potential 
Figures in brackets indicate tl terature references at the end of this paper measurements. 
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of the dropping-mercury electrode, the silver 
silver-chloride working electrode (Anode), and the 
sidearm of the saturated calomel reference electrode 
(Ref). The Ag-AgCl anode, of the thermal-electro- 
lytic type [3], had a large surface area to minimize 
polarization effects under the conditions of the 
measurements. This electrode was because 
it was found that the silver-wire anodes frequently 
recommended were of insufficient surface area to 
prevent polarization effects. The saturated calomel 
electrode (S. C. E.) was so constructed that a large 
portion of it was immersed in the bath. The side- 
arm contained an agar plug saturated with KCI 
All potential measurements were made with respect 


used 


to this reference electrode 

A number of capillaries were used during the course 
of the investigation. A capillary of marine barometer 
tubing producing drops having a mass of 6.9 mg 
mt value) was used in the investigation of the 
dependence of half-wave potential on concentration. 
In the determination of the dependence of the half- 
wave potential on capillary characteristics, the capil- 
laries had mf values of 6.8, 8.8, 10.2, 32, and 64 mg, 
respectively. The first three of these were of marine 
barometer tubing and the last two were of the con- 
Tavlor and Smith [4] 


stricted type described by 
the cell was cleaned in 


Before each series of runs, 
a hot sulfurie-acid—potassium dichromate solution 
and then steamed for about 15 min The stopper 
assembly was rinsed in distilled water and allowed 
to dry before it was placed in the cell The Ag-Ag( ‘| 
electrode was stored between runs in 0.01 .V HC] 
and before use was rinsed with distilled water and 
dried w ith a filter paper The calomel electrode was 
stored with the sidearm immersed in saturated KCI 
and was treated in a similar manner before insertion 
in the cell 

Nitrogen 
solutions. 


was used to from the 
The gas was passed over hot copper to 
remove anv traces of oxvgen before introduction 
into the cell While the measurements were being 
made, the nitrogen flowed over the surface to pre- 


remove Ooxvgen 


vent oxvgen from reentering the solution 

All measurements were made in a constant 
perature bath maintaimed at 25+0.01° C 

Stock solutions of the metal ions used were made 
by dissolving weighed portions of the pure metals 
in acid, evaporating the solution to dryness several 
times to remove all of the acid, and then dissolving 
the residue in distilled water and diluting to the 
desired concentration. <All materials were 
reagent grade chemicals without further 
purification. 

The  capacitance-conductance 
shown in figure 3 was assembled to 
resistance of the cells and the change of resistance 
with applied voltage. drop time, ete R, and R 
are the ratio arms of the bridge and were equal to 
500 ohms for the sensitivity used. The measuring 
arm of the bridge consisted of R, and C, a mica 
condenser variable in 0.001-uf steps from 0 to 1.11 uf. 

A 1,000-ohm resistance, Ry, was inserted to reduce 
the shunting effect of the potentiometer (P) on the 
detector (D). (P) was the slide wire of the model 


tem- 


other 
used 


bridge network 
measure the 
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Figure 3 Circuit for the measurement of the resistance of th, 
polarogra phi cel 


XI polarograph with the voltage increased to 5 
v. The signal source was a beat frequency oscillatoy 
with an output voltage hecessary to vive 0.04 
ncross the cell, The voltmeter (Y used to measure 
the d-e voltage drop across the cell was a Beckman 
H 2 pu meter. The detector 1)) was a cathode- 
ray oscillograph connected to the eireuit so that it 
could be used as an a-e voltmeter across the cell. 
as well as the null point detector for the bridge 
To study the variation of cell resistance with time. 
the bridge was set at some value greater than the 
minimum resistance and the time was 
from the beginning of the drop until a balance was 
indicated, This was repeated at oa series of R 
values until the curve shown in figure 4 was obtained 

Measurements of the half-wave potential, £ 
were made in the following manner. The current 
and voltage were measured at a point about 0.1 4 
or more preceding the half-wave potential to estab- 
lish a zero point for the wave, and at appropriat 
voltage increments up to and including the diffusion 
Measurements on the supporting 
corresponding voltages CrVE 
current, which were sub- 
tracted from the current readings above to. give 
values of ; Kor values sey ral tenths of a volt 
greater than the half-wave potential, the current 
is practically independent of the applied voltage 
(diffusion current plateau) and is designated the 
diffusion current, 7,. From these values of the 
current, log 4 (j¢—7) was caleulated and plotted 
with respect. to the measured voltage From. this 
plot, the value of / at log 7 (4,—7)=—0, when cor- 
rected for the voltage drop in the cell, corresponds 
to the half-wave potential. A typical plot is shown 
in figure 5. It may be seen that the curve ts linear 
and that the reciprocal of the slope is 0.0306. v, in 
satisfactory agreement with the theoretical value of 
0.0295 v. The value of R used in the correction 


recorded 


current plateau 
electrolyte alone at 
values of the residual 


for the 7R drop is four-thirds of the minimum cell | 


resistance corresponding to the average FA curing 
the life of the mercury drop. 
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a er 3. Results 
) 420 | 
= | 3.1. Variation of E;, With the Characteristics of the 
a Galvanometer Used on the Potentiometer 
v | 380 | 
— Measured values of the potential of the dropping 
electrode with respect) to a reference electrode 
_ { | vary with time during the growth of the mercury 
— drop due to the variation of voltage drop (/R) in the 
$ \ cell It is not feasible to measure either the maxi- 
5 mum or minimum value of this potential. Accord- 
- ingly, the potentiometer is balanced at some average 
ay c value as evidenced by the symmetrical deflection 
T 2 of the “null” galvanometer around its rest point 
: It is evident that such a potentiometric reading 
3 will be influenced by the characteristics of the 
— ’ null galvanometer and will approach the time average 
of the value when the characteristics are such that it can 
integrate the voltage variation that exists. 
It was found that for the conditions given below, 
when the damping Was near critical and the period of 
to & the galvanometer was of the same order as or greater 
“ae than the drop time, there was no significant variation 
4 \ of k with period ol the valvanometer or with the : 
nol - damping. If the period was short or the damping « 
‘man a ° ; . ‘ 6 resistance Was very large, the values obtained for the 
adie, TIME . SECONDS half-wave potential were increased. | nder these ° 
at it 2-05 0.6 57 5.8 0.9 0 conditions the galvanometer does not integrate the 
1] TIME ,(SECONDS)~'/3 voltage-time curve [5, 6], and the voltage read is no . 
ree we Sisctrical resistance of a drumuinarmerdery aletveds longer the average voltage corresponding to the . 
ee : oan . : | Pes Agate ag average current. The conditions for the measure- . 
: a 4 co eee ments were as follows: Drop time, 5.7 see; reducible . 
rded b) Det STasnO oe ‘ici sai ion, | mM Cd; supporting electrolyte, O.1 .N He I, . 
yer 0.01 pereent gel; =25° C; head 70 em; marine . 
if R barometer capillary ; reference, ae C. I. The val- 
ined vanometers and characteristics are listed below. 2 
ID j 
rrent ’ 
1 p Critical 
stab- ae NBS Rit Dested damping Damping resistane 
riate Vanomete! resist- 
ISlOn ance . 
rting | f ; 
AVE sec Ohms Ohms 
sub- 3008 L So 43 10 435 
wiv d 16350 3 83 2,700 700; 2,700: 
a 10442 5. 26 13, 000 10,000; 43.000 
volt . 5705 15.5 21,000 = 21,000; 
rrent 
ltage ? 
the | 
the 4 The magnitude of the error from this source ts not 
ited large but is significant in precise measurements. The 
this maximum error to be expected is of the order of one- 
cor- a f half of the 7? correction and would be obtained when 
onds using a balance indicator of very short period, for 
1own example, an oscilloscope. When using the galvanom- 
neal : eter with period of 1.85 see and infinite damping 
v, m p resistance, errors in the half-wave potential of 1 my 
a | ; ; ge were observed. On the other hand, with the galva- 
1 , mv nometer of 15.5-see period, half-wave potential 
sins Veevan & Piet of potential of the DME vevene ten Gai measurements were practically independent of the 


from which the half-wave potential is obtained, damping resistance used. 
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3.2. Measurement of Resistance 


found that, as a practical matter, the 
resistance measurements needed for the 7? drop 
correction may be made with sufficient accuracy on 
the supporting electrolyte, by the use of a small 
1,000-cps conductance bridge with an “electronic- 
eve” detector. provided that the Capacity balance 
maintaimed As shown in figure 4, the 
minus one-third power of 
The average value of 


lt Was 


is closely 
resistance varies as the 
time, in accord with theory. 
the resistance is four-thirds of the minimum value, 
occurring at the time of dislodgment of the drop. 
Measurements of with no applied d-e 
potential were found ‘o be applicable over the range 
of potential used, since the variation of resistance 


resistances 


with potential is small, and the accuracy of measure- 
ment decreases with the potential The resistance 
was found not to vary significantly with drop time 
For example, with a capillary for which m/—8.S at 
3.92, 3.12, 2.58, 2.22, 1.94, and 1.74 


drop times of 
196, 197, and 


sec, the resistances were 195, 196. 197 
197 ohms, respectively 

With i reducible ion in the solution, the behavior 
of thi 
potential wis venerally in accord with the observa 
However 


resistance and capacitance at the half-wave 
tions of other investigators |7. S as this 
phenomenon had ho direct bye aring on this investt- 
gation, it was not pursued in detail 


3.3. Variation of E With the Concentration of 
Reducible Ion 


Measurements of half-wave this 
study were made with a marine barometer capillary 
with a drop weight of 6.9 mg and at a drop time of 
The temperature was held at 25 —0.01° © 
for all runs. Measurements made on Cd in 1.0 
V ACI and 0.1 N HCI, and on Pb in OLN KCI 


and 0.1 N HCl show no significant variation of - 


potentials lol 


5S see 








with concentration of reducible ion from 0.5 mM to 
5.0 mM However, the correction for the 7R drop is 
very important in these measurements beimeg as 
large as 5 mv for the most concentrated solutions 
A typical plot in figure 6 shows curves before 

394 | 
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FIGuRE 6 Variation of E vith concentration | ind the 
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A) and after (B) the 7PR correction has bee: inade. 
lt was found that moving the reference e| trod 
from one position to another in the cell did not affee 


the potential measured with respect to the dropping. 
mercury electrode. A separate experiment was per- 
formed, using a cell in which three compart ments 
filled with 0.1 NV KCL were separated by tubes cop. 
taining agar plugs of several hundred olims resist. 
ance. Silver-silver-chloride electrodes were placed 
in the end compartments ana 
plied. It was found that the portion of the fall of 
potential between the electrodes, as measured by Q 
twee! 


2 potential Wiis ip- 


reference ¢ lectrode, depended on its distance be 
them In this case, however, the resistance at the 
surface of the Ag-AgCl electrodes is small, 
that of the intervening medium is large: therefor 
the greater part of the potential drop occurs across 
the solution and not at the electrode interface I) 
the case of the DALE, the resistance of the drop 
interface is large and that of the intervening medium 
greater part of the 
It is not 
possible, therefore, lo place the reference electrode 
close enough to the drop to exclude measuring this 
fall of potential Kou this rendsoth the half waive po- 
tential must be ih drop, althoug! ho 
cell current is involved 


Whereas 


is relatively small, so that the 


potential drop occurs across the mtertace. 


corrected for 


3.4. Variation of E:, With Capillary Characteristics 


reduction of a metal ion to metal solubli 


In the 


In mercury, the 
to vary with the drop time and size of drop of th 


hound 


half-wave pote ntial has been 








electrode, other factors remaining constant his 
Variation is approximately linear with drop time for 
a given electrode is is shown 1th fivure 7 lLlowever 
| ,’ 
p) 
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| 6 s 
4 } 
= 
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Figure 7. J on of I vith cap j 
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TABLE 2 
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of homogeneous electrode reactions in- 
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between 
several electrodes studied. 
have 
this investigation with capillaries having mf 
$2, and 64 mg, respectively. 


the 


been 


Half-wave pote ntial of Pb for a series of capularve s 


rhe supporting electrolyte is 0.1 N KCl] and 0.01 percent gelatin, except where 


oted. (Temperature, 25.0 











Capillary (mt value Drop time Fs 
mg sec t 

a8 1. 92 0. 3882 

&.8 2. 00 3891 

S.5 4.30 3906 

x8 5.74 $923 

32 3. 20 3918 

32 4.32 3917 

32 4. 86 3921 

32 6. 38 3927 

64 2.28 3913 

4 3.42 3922 

4£.9(0.1N HCI 8 3953 

6.9 5 SS8S5 
» 7254 " 
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d ee | 
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Ficure 8. Variation of Ey. with capillary characteristics in 


the case of the 
percent gel. 
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Fe > system ain 
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volving reduction from a higher to a lower oxidacion 
state, in which the dropping electrode acts as an 
indicator electrode, there is no variation with drop 
time and only slight variation with drop size, as is 


shown 1 


n figure 8. 


The supporting electrolytes were 


for Pb, Cd, and Tl, 0.1 NV KCI, 0.01 percent gel; for 


ke. 1.0 
0.04 


2, 3, 4, 


TABLE 3 


Phe suppor 


MZ Na, HPO,, 


values obtained for all cases are shown in tables 


M KCO, 


and 5 


Hla P-ane pote niial ol Tl lo 
and drop limes 


ting electrolyte 


» 


Capillary | 
nt value 


wo. 2 
10.2 
10. 2 


$2 
$2 
$2 


0.027 


uv ¢ 


p> Cirnie 


Se ¢ 


28S 
4.82 
». 4 


2.14 
$1 


5. 34 


a ser 


s0.1.N KCland 0.01 percent gelatin 


Trad 
Whew 


4519 
ho 
4548 
1548 
4584 
4587 
4593 


HH) 


4593 
4508 


ol 


, 0.0L percent gel; for quinhydrone 


Zo Nall, PO, 


l 


Half-wave potentials of the Fet*+, Fet++ system 
for a series of capillaries and drop times 


TABLE 4. 


he supporting electrolyte is 1.0 M K2C204 and 0.01 percent gelatin 
(Temperature 25,0° C 


Capillary . 
. ) . ; 
(mt value) Drop time Ess 
mg sec v 
10.2 2.14 0. 2321 
10.2 3. OS 2328 
10.2 4.26 2320 
32 3. 28 2333 
5. 42 2331 
64 2.10 2333 
4 3. 92 2333 | 
| 
| 


Half-wave potential of the Quinhydrone system for 
a series of capillaries and drop times 


TABLE 5. 


rhe supporting electrolyte is 0.04 M NasHPO, and 0.027 M NaHaPO, pH=7.0, 
(Temperature 25.0° C 


Capillary . 
3 ) . 7 

(mt value) Drop time Ei, 
mg sec t 
10. 2 2. 06 +0). O510 
10.2 2. US +. O518 
10.2 4.18 +. O507 
2 3. 24 +. 0408 
2 5. 38 +. 0493 
4 1. OS +. O500 
4 2. 0S +. 049090 








4. Discussion 


The results of this investigation show that the 
half-wave potential may vary significantly with the 
characteristics of the dropping electrode. Because 
variations of several millivolts can occur, it is impor- 
tant that the characteristics of the electrode be 
reported, in addition to information concerning the 
supporting electrolvte and temperature, to give the 
measurements their full significance. The impor- 
tance of the correction for the 7/? drop in the cell is 
also emphasized 

Strehlow and Stackelberg [9] have predicted vari- 
ation in /), with drop time caused by stirring effects 
within the Hg drop and the resultant disturbance of 
diffusion conditions for the amalgam. The variation 
of /. observed in this investigation in the case of 
the reduction of metal ions to an amalgam, anad the 
lack of such variation in oxidation-reduction svstems 
are ip accord with this theory The order of macni- 
tude is such that a few pereent change in the con- 
centration of the amalgam at the surface would 
account for the variation observed The problem 
of a quantitative treatment of this phenomenon Is 
under study 

lt has been shown bry other investigators | 10] that 
the half-wave potential is dependent on the capillary 
of irreversible reactions. In the present 
versus 


in the case 
case, however, the linearity of the log 7G P 
I. plot, and the pear theoretical slope obtained, 
attest to the reversibility of the reactions studied and 
eliminate irreversibility as the cause of the variations. 

In order to attam reproducibility within +0.2 my 
in half-wave potential measurements, several sources 
of error must The temperature of 
the cell and of the reference electrode must be known 
and held constant to within 01° €C 
temperature coefficients are approximately O.S my °C 
in each ease. The residual current must be deter- 
mined from measurements on a cell containing sup- 
porting electrolvte only, and the corresponding 
correction applied to the limiting current at each 
point. Fortunately, moderate errors in the residual 
current can be tolerated, as it is only a small fraction 
of the total current in the usual case 


be considered. 


since the 


The large st sources of error arise in the ev aluation 
of the versus E relationship and in the 
determination of the voltage drop in the cell With 
respect to the first. comparison of graphical methods 
with least squares adjustment of the data shows 
agreement in the two methods to 0.1 mv. The 
simple plotting procedure ean thus be reled on 
within this limit of precision. The magnitude of 
the 72? correction depends on both the current and 
the cell In measurements of half-wave 
potential, cells should be designed to offer minimum 
electrical resistance, and the design used in this work 
appears to approach the optimum from this point of 
Such cells have resistances of the order of a 
few hundred ohms (largely interfacial 
which can be readily measured with ordinary equip- 
The magnitude of 


log } } } 


resistance. 


view, 


resistance 


ment to within 5 to 4 percent 
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the diffusion current depends on both the siz 


i Of the 
dropping electrode and the concentration if the 


electrolyzable material. However, since the resist. 
ance of the cell decreases with increase of the size of 
the drop, with corresponding increase in the currepy 


it is not difficuli to choose a combination so t! ut the 
uncertainty in the 7? correction will not exceed +0) 
mv, except in the case of the largest electrodes, wher 
errors of several times this amount may result 

An Important source of error in the measurement 
of half-wave potentials is the junction potential of 
the reference electrode-solution interface 
may be of the order of several tenths of a millivol; 
and has been recently discussed by Vleek {11 No 
corrections for junction potentials have been Tack 
in this study, as they would be only approximat 
and would not affeet the precision, 


This erro; 


5. Summary 


The sources of error in half-wave potentials have 
been investigated and found to he chiefly in the cor- 
rection for the 7? drop in the solution and in the 
actual reproducibility of current and voltage meas- 
urements. When proper precautions are exercised 
it is possible to obtain values reproducible to within 

0.2 my 

The half-wave pote ntinl is shown to varv with the 
size of drop and with the drop time of the electrode 
in the case of the reduction of metallic tons to metal 
soluble in the mereury electrode On the other 
hand, much less variation was observed in homoge- 
neous electrode reactions ino which the dlropping- 
mereury electrode 
Accordingly, the observed variations are attributed 
to stirring effects within the drop and the resulting 
disturbance of diffusion conditions for the amalgam 
formed It) heterogeneous rencriions., In View of thi 
large variations found, values of half-wave potential 


acts aus an indieator cleetrode 


are not significant unless the characteristics of th 


electrode used for the measurements are reported, 
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Infrared Measurements 


From 50 to 125 Microns_ 


Earle K. Plyler and Nicolo Acquista 


\ prism spectrometer has been changed 
removing the prism and by replacing the Li 
and ISO lines per inch were used for 
ervstal quartz window served as the detector 


isc of reflection filters Three types of re fle 
were gratings with 2,000 lines per inch, rot 
of NaCl, KBr, and CsBr The absorption 
the spectrometer Phe spectra of several 

Ineasured, and absorption bands have been f 


1. Introduction 


During the last 2 vears prisms of Csl have been 
available, and it is now possible to make measure 
ments 50 gw [1] (As quartz has its strongesi 
reststrahlen band at 9 uw, it should) transmit: radia- 
tionat LOO w. For thicknesses less than | mm, quartz 
has a good transmission of infrared radiation im the 
However, with prisms of 
was 


lo 


long wavelength region 
erystal quartz of a refracting angle of 7.5°, it 
found that no appreciable energy was transmitted in 
the LOO py Mnerg, was recorded in the 
150-u region, but the resolution was low. Diamond 
is another material which is transparent in the long- 
wavelength region, but the expense and difficulty in 
obtaining material of sufficient size for a prism for 
infrared measurements would make its use impracti- 
llowever, windows of diamond for detectors 


region OU to 


cable 
are useful for the infrared region 


2. Experimental Procedure 


In order to extend measurements bevond Ow, 
experiments were carried out by the use of small 
gratings.  Large-grating spectrometers for the far 
infrared have built in several laboratories [2], 
and it was felt that sufficient energy would be avail- 
able with small gratings. A graiing was installed 
ina spectrometer in place of the Littrow mirror with 
the prism removed. At first the observations were 
made with an 1S° off-axis collimating mirror, which 
was rotated a few degrees in order to fill the grating 
with the radiation. A small improvement in the 
resolution was obtained later on by installing a 14 
Of-axis mirror. ‘Two gratings were used in the 
spectrometer, which had ruled surfaces of about 24; 
by S inches and were ruled with 320 lines per inch 
and 180 lines per inch. The spectrometer is able to 
accommodate slightly larger gratings, whieh would 
give an increase in energy and resolution. The 
grating was placed in a holder which allowed ad- 
justments of its position, so that a mercury line in a 
higher order could be brought sharply in focus at 
the exit slit. , 


i work wa lone under 4 procram of basie instrumentation, ipported n 
rt by the Office of Naval Research, Air Research and Development Command, 
ithe Atomic Energy Commission 

Figur brackets indicate the literature references at end of this paper 


the region 50 to 


ong wavelength grating spectrometer by 
grating Ciratings of 320 

\ thermocouple with a 

Stray radiation was greatly reduced by the 
‘ction filters were found to be effieient They 
ighened silver mirrors, and reststrahlen plates 
vapor were used to calibrate 
ethvlenes and ethanes have 
50 and 100 microns 


tow 
ttrow mirror witha 
125 microns 


lines of water 
substituted 
ound between 


bee if 


The curvature of the entrance slit was not propet 
for the grating, but this did not reduce the resolution 
much because slits of 0.5 mm or larger were used for 
the region from 50 to 125 4 Hlowever, better resolu- 
tion will be obtained with straight slits. A thermo- 
couple with a Csl window was used from 2; 
and another thermocouple with a crystal quartz 
window was found to be suitable from 50 to 125 4g 
It is planned to use a thermocouple with a diamond 
window from 25 to 125 yu, thus avoiding the changing 
of the thermocouple for the different regions of the 
spectrum 

Special attention is necessary to the problem of 
stray radiation in the long-wavelength region. The 
stray radiation was reduced to a few percent by the 
use of reflection filters and = reststrahlen plates im 
place of the mirrors in’ the spectrometer. The 
optical svstem of the grating spectrometer is shown 


9 to 56 pw, 


in figure | after the grating had been installed and 
the plane mirrors replaced by reflection filters. The 
two small plane mirrors, immediately in front of and 


behind the exit slit: of the spectrometer, were re- 
placed by two plane gratings, and the plane mirror 
in the source box of the spectrometer was replaced 
by a reststrahlen plate of NaCl, KBr, or CsBr 
Also, three gratings could be used in place of the 
mirrors, but the reststrahlen plates are very useful 
in reducing second order spectra. A system which 
was as effective as the gratings and the reststrahlen 
plates in the removal of stray radiation consisted in 
the use of blackened silver mirrors which were first 
roughened with fine emory cloth. These reflection 
filters have been deseribed previously [1]. 


Figure | The oplical system of the grating spectrometer, 


A, B, and C are reflection filters which have replaced the mirrors of the spec 
trometer. G is the dispersion grating which is placed on the Littrow mirror 
mount, and ‘T is the thermocouple with a quartz window 
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\ study has been made of the efficiency of different 
types of filters for removal of stray radiation. The 
water-vapor spectrum from 45 to 70 uw was observed 
by the alternate use of several reflectors for the 
filter A of figure 1, the other parts of the 
trometer being unchanged during the comparisons. 
Figure 2 shows the results obtained for two gratings 
with 1,000 lines per inch, one roughened silver mirror, 
and a plate of NaCl. The zero line represents the 
recorder reading when an opaque shutter ts placed in 
the beam. The water-vapor absorption at 139 
em~', under the atmospheric conditions of the ob- 


spec- 











servations, was 100 percent. By measuring hoy 
much the bottom of this absorption line is above thy 
zero line an estimate of the stray radiation can }y 
obtaimed. An estimate of the amount of the desired 
radiant energy reflected from the filter can be deter. 
mined by measuring the difference between the de. 
flection at the minimum of the line at 139 em~! and 
the maximum at 163 em Grating 1 is the most 
efficient of the four filters tested. It has 1,000 lines 
per inch and the plane surface was almost completely 
cut away in the ruling so that the central image had 
about | percent reflection. The other vrating 





Ficure 2. 


The base line for each filter represent 
I 


NUMBER 


A comparison of the efficiencies of different types of reflection filt rs. 


the deflection with an opaque shutter in the beam 
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gratiig 11) had 1,000 lines per inch and the reflection 
for normal ineidence was about 5 percent. Grating 
11 is much less efficient in reducing the stray radia- 
tion Phe roughened silver mirror does not reduce 


the desired energy appreciably but does allow con- 
siderable stray radiation to enter the spectrometer. 


The stray radiation can be made less for the silver 
filter by producing a greater amount of roughening on 


the surtace, but this also reduces the desired energy. 
Reststrahlen plates of NaCl, KCl, KBr, and CsBr 
are eflicient filters for the region from 45 to 125 u. 
The spectral region of high reflection is limited to 
narrow bands and one material will not cover a wide 
range of the infrared spectrum. The CsBr 
strahlen plate is useful in measurements from 100 
to 125 u. The second-order spectrum is intense in 
this region from the 180-lines-per-inch grating and it 
is not removed by the grating filters or the roughened 
silver mirror. 

Since the quartz window on the thermocouple ab- 
sorbed the radiation from 5 to 40 uw, the second order 
from the grating did not overlap the first order 
except in the region from 90 to 125 uw when the 
roughened silver reflectors were used. 

With these changes of the mirrors in the optical 
system just mentioned, it was found that sufficient 


rest- 


energy for measurement could be observed from the 
globar source with the same amplifier that was used 
for the near infrared region. 

The instrument has been tested and calibrated by 
the use of the water-vapor spectrum in the region 
from 25 to 125 4. It has been found that the instru- 
ment can be used in three ways. First, by the use 
of the second and third orders of the 320- and 180- 
lines-per-inch gratings, high resolution can be ob- 
tained in the region from 2: Second, the 
320-lines-per-inch grating, with the thermocouple 
having a quartz window, gives good resolution in 
the first order from 50 to 70 uw. Third, the 180- 
lines-per-inch grating can be used from 50 to 125 u. 
Gratings with various spacings can be used to obtain 
the best resolution in a particular region. One grat- 
ing with 240 lines per inch would be suitable for the 
region from 50 to 100 uw. 


5 to 50 M. 


3. Results 


Figure 3 shows the water-vapor absorption lines 
from 23.5 to 33.5 uw as observed with the 320-lines- 
per-inch grating. The top curve is the second-order 
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Figure 3. The water vapor spectrum recorded with a 320-lines-per-inch grating. 


The upper curve is the second-order spectrum and the lower curve is the third-order spectrum. 


occurs at 305 em='. 
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and the bottom curve is the third-order 
spectrum. ‘The five lines in the region of 395 em 
are well resolved by the instrument. The strong 
absorption line at 340 em”! was not reported in the 
earlier experimental work for this region, and this 
transition is predicted from the energy levels in- 
vdlved to have appreciable intensity. In the region 
of 300 cm™! the third-order spectrum is oveslapped 
by the second-order spectrum; this overlapping can 
be avoided by the use of a NaF reflector for this 
region 

Figure 4, upper curve, shows the region from 50 to 
65 was observed with the 320-lines-per-inch grating, 
and the lower curve is for the region from 65 to 100 u 


spect! 1}, 


which was measured with the 180-lines-per-inch 
grating Also the region from 100 to 125 uw. measured 
by the use of a CsBr reflector, is shown in figure 5. 


Slits of 1.8 mm were necessary for the measurements 
in this region . 

The lines drawn above or below the spectrum 
represent the positions of the water-vapor lines as 


calculated from their energy levels These wave- 
numbers are reduced to vacuum and have been 
published by Benedict, Claassen, and Shaw [3]. 


Only single lines have been used for the 
calibration, but the approximate positions of some 
complex lines are also listed on the figure. The lines 
were identified by using the wavenumbers obtained 
from the grating constant. While the Wwater-vapor 
spectrum is Of great value for calibration and testing 
of the imstrument, it should be removed for the 
measurement of other spectra. This can be best 
accomplished by placing the spectrometer in a vac- 
uum chamber. The absorption of the water vapor 
was ereathy reduced hy circulating dry nitrogen in 
the enclosure of the spectrometer, In the measure- 
ment of the absorption bands of liquids and vapors 
the removal of the water-vapor spectrum is essential 
to obtain band shapes 

The primary purpose of developing this instrument 
was to have available a spectrometer for measuring 
the absorption spectra of materials in the region 
from 50 to 125 wu. In the study of the vibrational 
bands of many molecules it was found that it was 
difficult to make a complete assignment of the differ- 
ent modes of vibration because several of the bend- 
ing vibrations would be between 50° and 
Since the instrument has been in operation the long- 


strong 


125 yp. 


wavelength spectra of C,ClL,, C.Bry, CLF, Br have been 
measured and their spectra analyzed. Cells with 
windows of erystal quartz, 0.2 mm in thickness, 
were used for these measurements. Many of the 
spectra of the substituted ethylenes, ethanes, and 
benzenes which contain fluorine, chlorine, bromine, 
or iodine will contain bands bevond 50 uw, and an 
instrument of this tvpe will be of great value in 
obtaining sufficient data for a complete vibrational 
analysis. As an example of measurements on long- 
wavelength spectra figure 6 shows the absorption 
bands of CH.CI-CClL,, CHBr,.-CH,Br, and C,HC! 
between 50 and 100 4. The longest wavelength 
band which was observed for these materials aup- 
peared at 114 em™!. With a cell of 1 mm thickness 
it was necessary to dilute CH,CI-CCl, with carbon 
disulfide in ratio | part to 6 so that the bands could 
be resolved 

The results that have been obtained demonstrate 
that spectra can be measured in the region from 50 to 
125 uw with small gratings. For the arrangement with 
one grating of 200 to 240 lines per inch and a thermo- 
couple with a quartz window, the costs should not 
be greater than $1,000. Where several gratings are 
emploved and a vacuum enclosure with the necessary 
external controls is made, the cost would be greater 
The number of gratings which will be employed and 
the changes in the filters for different regions will of 
necessity depend on the nature of the measurements 
which are being made. For the measurement of 
Vibrational spectra and the transmittances of solids 
and liquids a single grating is adequate. 
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Surface Tension of Molten Glasses in the System 
Barium Oxide-Boric Oxide-Silica 


Herman F. Shermer 


BaQ-BoO 


irface tension of 154 compositions in the system -SiQ, were measured 
bv a maximum pull-on-eylinder method, over the approximate range 900° to 1,300° ¢ 
The region of glass formation was investigated along with the adjacent regions Several 
compositions were studied in the two-immiscible-liquid region. In this two-liquid region 


the lines of constant surface 


the two-liquid region the surface tension was 
lesser degree by SiO An 900" © 
corresponding to the 3Ba0-3B.0 


1. Introduction 


A study ol the surface tension of the system 
BaO-B.O.-S10, was undertaken as part of a compre- 


1 


hensive investigation of various physical properties 
svstems |] to 135}. The 
system, BaO-B,O,, has 
been reported of a previous investigation of 
the properties of the binary alkaline-earth borates 
10) The studied 
limited by relatively bigh 
region of high silica. The phase 
the svstems BaO-B.O, and BaOQ-P.O,- 
FT| two-liquid region at low concentrations 
of alkaline-earth oxide 
The svstem BaOQ-B.O fundamental 
importance in the ceramic field and especially in the 
starting poimt for 
Class of berium-crown 


of simple glass-forming 
surface tension of the binary 
is part 
range ol compositions 
ervstallization and by 


Wils 


yiscosityVv it the 
diagrams [ol 
SiM, show 


si) is ol 


class industry, as it 


serves 2s a 
investigations of a large 


glasses 
2. Experimental Procedure 
2.1. Preparation of Glasses 
The composition of the glasses used in this investi- 


shown in table J 
13oo lo 


The glasses were made 
1.500 (" the required 
chemicals ino platinum erucibles and 
stirring the melts with platinum stirrers until they 
appeared homogeneous 
Irom 100 to 500 2 
Chemical-analysis samples were taken from the 
melts before and after surface-tension measurement 
The reported compositions are the average of the 
two determinations unless otherwise noted in table | 
The BaO and SiO 
metrically 


fallon is 
by melting at 


amounts ol 


The melts varied in weight 


contents were determined CTravi- 
wand the BLO, by difference. Small 
changes were noted between the batch composition 
and the first analysis \ further loss of a few tenths 
percent of BLO, was noted between the two 


of i 


tension followed the tie 


2510) compoul d, which 


the phase diagram Outside 


addition of BaO and to a 


lines in 
increased by the 


the surface tension was relatively high near the composition 


has a melting point of 1,009° C 


analyses, indicating volatilization during the tests. 
The estimated probable errors of the gravimetric 
0.1 weight percent, 


2.2. Method of Test 


determinations are about 


The method of measuring surface tension used in 
this investigation has been previously described in 
detail [14]. It consists in determining the pull ex- 
erted by the test surface on a thin-walled platinum 
evlinder during the removal of its lower edge from 
contact with the liquid. The first determination of 
the surface tension of any melt was made at 1,300° C 
Subsequent measurements were made at 100-deg C 
intervals, decreasing the temperature until the glass 
Was too viscous to yield in a reasonable time or until 
cry stallization occurred 

The density values used in the correction of the 
calculations were determined by 
interpolation between room-temperature values for 
and a few known values in the 
temperature and composition range of interest. For 
convenience in presenting the data, all values of sur- 
face tension have been adjusted to even 100-deg 
intervals. The temperature of any measurement 
may have been 4 deg C from that given in table 1. 


surface-tension 


each composition 


3. Results and Discussion 


Surface-tension data are given in table 1. Differ- 
ences in batch materials and source of glasses are 
indicated in the footnotes. Figure 1 shows the iso- 
epitatic 15], or constant surface-tension, lines at 
1,300° © and the locations of the compositions in- 
vestigated. The surface tension increased with de- 
crease in B.O, content. This figure also contains an 
outline of the two-immiscible-liquid region and the 
conodes, 1. e., the tie lines connecting the two con- 
jugate compositions in the phase diagram as deter- 
mined by Levin and Ugrinie [16]. 

In the two-liquid region the isoepitatic lines ap- 
proximate the conodes. This is in accord with the 


Che sampl were analyzed | Ralph I. Cohen, Ernest B. Clark, and J. 
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BaQ-B.O.-Si0O f J.3000> © and locations of compos 


ca 


fact that at any given temperature the upper laver 


i all CODLPOSLULOs along a conode Is identical, The 
ighter liquid will tend to rise to the surface, but if 


there were Incomplete separation, there may be some 


if the lowe! licytaicl in the upper laver The com- 


pleteness of separation becomes greater as one de- 


parts trom the two-liquid Interface: and at the sur- 


surface tension is measured, the 
minimum amount of the lower laver will be present 
The presence of lower liquid in the upper liquid is 
ndicated by opacitv. The glass that adhered to 
ne ring’ Was clear, indicating the probable absenes 
of anv of the lower liquid in the surface under test 
Therefore, surface measurements may give a good 
ndication of the direction of the conodes in the two- 
quid region, and surface-tension measurements may 
offer another approach to the study of phase equilib- 
ain the two-liquid region 

Figure 2 shows the isoepitatic lines at 1,100° © 
{ comparison with the previous figure shows a 
shiehit counterclockwise rotation of the lines within 
the two liquid region This is in aceord with the 
otation of the conodes with temperacure reported 

thre phase work 

Figure 3 shows a plot of the Isoepitatie lines at 
90° CC. The changes in surface tension between 
L100° and 900° Car 
L300° and 1.100° © 
The constant surface-tension lines near the com- 
3Ba0-3B.0,-2810., 
compound deviate from simple curves. This com- 
pound was characterized by Levin and Ugrinie 
l6] as one that melts congruently at 1,009° C and 
Whose primary field occupies a large portion of the 
glass-forming region of the system. This is 109 
deg C above the temperature of the present measure- 


similar to those between 
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| IGURE 3 I soe p latic lines of the system Bat -B.0O.-Si0O at 
900° C., 
he open « e indicates the location of the 3BaO0-3B 05-2810. compound 
tified by Levin and Uegrinie [16 The two closed circk ndicate the loca 
tions of the low-surface-tension areas 


ment. The deviation of the isoepitatic lines is associ- 
ated with an irregular surface-tension region near 
the compound (indicated by the open circle) and 
two relatively low surface-tension areas (closed 
circles 

Above the liquidus temperature (1,009° C) at the 
compound composition no departures of the isoepi- 
tatic lines were observed. The possibility of the 
formation of crystals may be the cause of these 
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4. Summary and Conclusions Jou 


The surface tension of a series of BaO-B.O.-Sid 
COMM positions Was measured over the range “OO to 
1.300° CC. In the two-liquid region the values 9 
surface tension and the coefficients of surface tension 
followed the conodes in the phase diagram, indir ating 
i possible new technique for phase study in this type 
of region. In general, the surface tension increase, 
with decrease in BL.O, content. At 900° C the lines | 
of constant surface tension reflected an irregularity 
near the composition corresponding to the 3Ba0 
$B.0,-2810, compound, which has a melting point 
0.020 of 1,009° C 
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Relationship Between Crystal Orientation and Stress- 
Corrosion Cracking in Alpha and Beta Brasses 


Hugh L. Logan 


The 


large-grained 


in large-grained beta brass were determined by X-ray diffraction methods 


crystallographic planes bounding intererystalline stress corrosion cracks in two 
Alpha brasses and those followed by transerystalline stress-corrosion cracks 


Intererystalline 


stress-corrosion cracking in the alpha brasses originated between grains that, because of 


their relative orientation, had high interfacial energy contents 


It is the relative orientation, 


not the erystallographie type, that determines the susceptibility to stress-corrosion cracking 


of alpha brass in grain boundaries approximately normal to the applied stress. 
occurred 
favorably for slip, and in planes that were approximately normal to the applied stress. 


talline stress-corrosion cracking of beta 


brass 


‘Transecrys- 
nost 


No 


in grains that were oriented 


plane or family of planes in beta brass was found to be particularly susceptible to stress- 


corrosion Cracking. 


mechanism 


1. Introduction 


Stress-corrosion cracking, or season cracking, in 
brass Was first reported in 1906 [7].! It has been the 
subject of many investigations and several compre- 
hensive bibliographies have been published [2, 3]. 
The study of the stress corrosion of brass is partic- 
ularly important from a mechanistic point of view 
because stress-corrosion cracking is predominantly 
intererystalline in the face-centered alpha brass and 
mainly transcrystalline in the body-centered, ordered 
beta brass 

If corrosion occurs in unstressed alpha brass, it is 
generally of the pitting type; only two reports [4, 4] 
of intercrystalline corrosion have come to the author’s 
attention. However, if the material is exposed under 
sufficiently high residual or applied tensile stress in 
an atmosphere containing ammonia and water vapor 
the specific corrodent for brass), predominantly 
intererystalline stress-corrosion cracks will develop. 

An excellent discussion of the mechanism of stress 
corrosion in alpha brass was presented by Read, 
Reed, and Rosenthal [6] at the ASTM—AIME Sym- 
posium on Stress Corrosion in 1944. Read et al 
proposed the hypothesis that the grain-boundary 
regions in this brass have higher energy contents 
than the grain interiors, and therefore the grain 
boundaries have a greater susceptibility to corrosion. 
These workers discounted the possibility of a grain- 
boundary precipitate but suggested that lattice im- 
perfections may explain the higher energy at the 
grain boundaries. In 1950 Robertson [7] suggested 
that relative grain orientation could be a factor in 
materials susceptible Lo intercrystalline stress-cor- 
rosion cracking. So faras the author knows, expla- 
nations of the mechanism of the transecrystalline 
cracking beta brass have not 
reached the same stage of development as those for 
alpha brass. 

Following World War IL a great deal of interest 
developed in the study of metal boundaries or inter- 
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It is postulated that in each case cracking progressed by a film rupture 


faces and the results of a number of investigations 
were summarized in a symposium, “Metal Inter- 
faces” [8], published in 1952. . In this symposium it 
was pointed out that interfacial or grain-boundary 
energy is a function of the difference in orientation 
between the grains [9]. A number of investigators 
have shown that the interfacial energy increased to 
maximum values with increased difference in angular 
orientation between pairs of crystals for angles rang- 
ing from 12° for tin [/0] to 29.8° for an iron-silicon 
alloy [7/7]. If the grain-boundary energy, and there- 
fore the chemical activity of grain boundaries in 
alpha brass, increases with increased difference in 
orientation of adjacent grains, it is to be expected 
that the probability that a stress-corrosion crack 
would develop at a particular grain boundary would 
depend on the relative orientation (or degree of 
misfit) of the crystals meeting at that boundary. 
Cracks should occur most frequently between crys- 
tals whose relative orientations were in the range in 
which maximum grain-boundary energies have been 
obtained for other metals and may be expected for 
brass. 

It is postulated that there is a greater probability 
of corrosive attack at the grain boundaries having 
the highest interfacial energies than at other points 
on the surface of the brass. The corrosive attack 
at these grain boundaries would form sharp notches. 
Those notches that were most favorably oriented 
with respect to the stress pattern in the metal would 
then act as points of stress concentration that would 
produce eventual rupturing of the protective film on 
the metal at the root of the notch. Stress-corrosion 
cracking would then proceed by a film-rupture 
mechanism, as discussed in an earlier paper [72]. 

As a part of the current study of the mechanism 
of stress-corrosion cracking, the crystallographic 
planes bounding stress-corrosion cracks at their 
origin and the angles between the crystallographic 
planes in the crystals on opposite sides of stress- 
corrosion cracks have been determined for specimens 
of two compositions of alpha brass. 
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It is generally accepted that, in the two-phase 
alloy, stress-corrosion cracks will be intererystalline 
in the alpha phase and transervstalline in the beta 
phase. The of brass has not 
been as extensively as that of alpha brass 
However, from the observations of stress-corrosion 
cracking in the two-phase system it has generally 
been considered that stress-corrosion cracking would 
be transerystalline in a single-phase beta alloy. As 
a part of a study to determine why stress-corrosion 
cracking intererystalline in’ alpha and 
transerystalline in beta the ervstallographie 
planes on which the transerystalline stress-corrosion 
cracks had their origin were determined for 
several specimens of large-grained beta brass 


corrosion beta 


studied 


stress 


was brass 


brass, 


also 


2. Materials and Experimental Procedure 


The two alpha brass alloys studied had the nomi- 
nal compositions of 90 percent C‘u-10 percent Zn 
bronze and 70 pereent C'u-30 percent 
Zn (cartridge brass). The beta brass had a compo- 
sition of approximately 50 percent Cu-50 percent Zn 
Ingots of these alloys approximately l by I by Sin. 
molds, had grain diameters ranging 
Sections containing equianxial erys- 
ingots for this investt- 
three alloy 5, 


commercial 


sand 
to |, in. 
tals were selected from these 
gation. In addition these 
large-grained specimens of the 90) percent Cu-10 
percent Zn alloy were prepared I> the strain-anneal 
method and were also used in this investigation 
Specimens having the shape and dimensions shown 


east in 
from ‘'s 


lo cust 


in figure | were machined from the ingots and the 
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rod. The flat surfaces of the specimens Were re. 
peatedly polished (using 000 metallographic polishing 
paper) and etched until Laue back reflection X-ray 
patterns indicated that the surfaces were free frop 
cold-work. The end sections were then coated wit} 
a waterproof lacquer, leaving only the central 

in. exposed to the corrosive medium. Th speci. 
mens were loaded in tension under a stress of 5.000 
to 6,000 Ib/in? and exposed in an atmosphere cop. 
taining ammonia, water vapor, carbon dioxide, and 
air until stress-corrosion cracks appeared. ‘Two gyp. 
faces, at right angles to each other, specimen of 
alpha brass containing stress-corrosion cracks are 


shown in figure 2. 


of ; 


To determine the erystallographic plane or planes 
bounding or followed by stress-corrosion cracks, the 
cracks studied should extend continuously from on, 
surface of the specimen into an adjoining surface ay 
right angles to the first. It is also necessary that thy 
two crystals bounding the crack be continuous, sg» 
that surfaces Ay and figure 3, be surfaces of the 
same crystal machined at right angles to each other 
Laue back reflection X-ray photographs were mad 
of the erystal surfaces A, and A, of crystal A and 
projections of these surfaces, relative to their 
respective cubic axes, were plotted on stereographi 
by Greninger [13] 
Was rotated through 


nets, using the method deseribed 
The plot of surface A,, figure 3, 
90 deg about the N—-S axis of the net, and the net was 


then superimposed on that for surface A If the 
cubic poles of the two nets coincided (within the 
experimental error), the surfaces A, and A, wer 
5/16"—I8 NC 5/i6" 
nw) 3 
a ULL) t 
oe— 7/8"—a— 3/4 5/16" + 
Z 
ve | 
f i 
FIGURE 3 Schematic drawing of « pstala A and B38 with stress- 
corrosion crack indicated hy arrow mak nq an angle @ wil 


corner of specimen in face l and angie @in tace 2 
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eonsice od to be faces Thhine hined normal to each other 
yvstal A. In the study of intererystalline 


from ¢ 
a similar check was made on crystal B 


erackin 
The an 
made with the common edge of the specimen were 
determined on surfaces | and 2, respectively, with a 
These angles then being 
the traces of the stress-corrosion crack could 


@ and @ that the stress-corrosion cracks 


les, 


toolmak « r’s Microscope. 
knowl ; 
pe plotted on the stereographic projections for the 
wo surfaces of the crystals [74], and the normal to 
the plane of fracture could also be plotted on each 
projection. ; 

In order to determine the difference im orientation 
between the two crystals bounding intererystalline 
stress-COrrosion cracks, the angles between the cubic 
poles *and the normals to the planes of fracture in the 
wo crystals were measured using a Wolf net. From 
these data the normals to the plane of fracture in any 
erystal could be plotted in a single 100, 110, 111 
stereographic triangle, and if the fracture was inter 
crystalline the angle between the normals to the 
planes bounding the fracture was determined. 


3. Results and Discussion 


The corners of the specimens were points of stress 
concentration and examinations of the specimens 
indicated that stress-corrosion cracks originated 
there. ‘This facet and the rectangular cross section of 
the specimens made possible the determination of 
the ervstallographic plane or planes followed by or 
hounding the Stress-COorrosion cracks 


3.1. Alpha Brass 


The orientations of the planes bounding 13 inter- 
crystalline stress-corrosion cracks in 70-percent-Cu 





-percent-Zn brass and IL cracks in 90-percent- 
Cu-10-percent-Zn brass) were determined. Data 
vere Obtained from four different specimens of 
each alloy, each of which contained more than one 
ils | tted on the projection w surface 2 ofery ils A or Bomust be 
ct dence with those i inti 1 tev thee ame rotation that will 
\ ibic poles into comeidence In the case of intercrystalline fracture 
¢ projection for surface 1 of et tul A uperimposed on that for surface 1 of 
tal B hat their N-S pol pre musty aligned with a fixed point on the 
ecime! vitudinal axis), coincide, the normals to the fractured surface will 
i milar results are also obtained for surface 2. It is thus seen threat 
eisno ambiguity as to the normal to the plane of the stress-corrosion crack 
Phe cubie poles of crystal A are designated orthogonal axes Na, Ya, Za, and 
gles between the normal and these poles are designated /’r, Py, and ?:, 
pectively Angles determined on face 1, designated 7’,,, ete., were checked by 
m of angles on face 2 designated 2? ote Asa further check use 
the expre lon 
I? $ Pyol4 i’ / + cos Py cos Pay) +A<l 
{thea ‘ ve ver than O.O10 it indieated t! 
cking the measurement of the angk An average value for 2’,, ete 
} } ned as follow 
/ . P,, cos I \ 

il Be ild t rmined and designated Q,, ete 
between the erystallographic plane m opposite ides of the 
be nined unalyt | providing 2’,, Py, Os, Gy, ete., are 

! ler re ! a tude of the angles for each erystal) from 
p ‘ 
PrP, cos Q Pt Y, +'cos Ff’, cos Y, 
il Th ok Unie ue u wt! ingle between two crystallographic 
smnot be obtained Phe angle determined either graphically or 
ica by the methods deseribed above is the smallest angle between the 
plat 
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erack. The results, plotted ina 
single stereographic triangle, are shown in figure 4 
The crystal boundaries and the consequent inter 
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vranular cracks do not necessarily follow the same 
eryvstallographie planes throughout their entire 
lengths. Hence at different points on the boundary 
between two erystals the bounding planes may be 
different In this investigation the bounding crystal 
lographic planes were determined not only at the 
incidence of cracking but also at some other pomnts 
along the continuous path that was representative 
of the general path or trend of the crack between the 
same two erystals. Poles of the latter planes were 
distributed in the same general way as those shown 
in figure 4. The listribution of the normals to all 
of these planes both at the incidence and along the 
veneral path of the crack is random. These data 
definitely indicate that it is the angle at which the 
crystallographic planes intersect that determines the 
susceptibility of a particular grain boundary, in alpha 
brass, to the development of stress-corrosion cracks. 

The angles between the erystallographic planes in 
the paired crystals of the two alloys, determined by 
the methods deseribed above, are given in table | 
Values were rounded to the nearest degree and are 
given in order of increasing magnitude, 

Kor more than 70 percent of the cracks studied, 
the angles between bounding crystallographic planes 
at the origin of cracking were in the range of 13° 
23°, inclusive; the minimum angle for all the pairs of 
crystals studied was 8 Angles of 13° to 23° are 
in the range ino which maximum grain-boundary 
energies have been reported for the face-centered 
metals silver [/4] and lead [7/46], as well as for tin [/0| 
and an iron-silicon alloy {/7/] 

No experimental data for brass have come to the 
author's attention; however, a value of 22.2 
been reported for copper [171 


lias 
The data indicate 
that the intererystalline stress-corrosion cracking in 
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alpha brass is associated with grain boundaries that 
may be expected to have high energy contents. 

It is suggested that in an environment containing 
the ammonium ion, galvanic cells will be set up 
with the high-energy grain boundaries as the anodes 
and the erain faces or other lower-energ) areas as 
the cathodes. The circuits are completed through 
the interior of the alloy. The metal will be dis- 





solved at the grain boundary anodes and the resy}. 
ing pits will have the appearance, in cross ection, 
of sharp notches with their apexes on the grajy 
boundaries. On the exposed surfaces of the Crains 
the ratio of anodic to cathodic areas will general}; 
be much larger than those of the grain boundary 
areas to their cathodes. Hence pits formed on the 
exposed surfaces of the grains will be much shallower 
and will have a more rounded cross section than thos 
formed at grain boundaries. In the absence of 
stress, pitting in either case will become more or leg 
self-limiting, because of polarizing effects, and inter. 
crystalline corrosion extending to any appreciable 
depth will be extremely rare. On the other hand 
sharp notches centered on grain boundaries and favor. 
ably oriented with respect to either residual Or ap- 
plied tensile stress patterns will act as regions of 
stress concentration. As the notches become suff. 
ciently deep or the stresses sufficiently high, the 
protective film on the brass at the roots of the 
notches will be ruptured, exposing metal that has 
an electrochemical-solution pot ntial approximately 
ly v more negative than the metal covered with a 
atmospherically formed film or with  corrosio; 
products [12]. The unfilmed, or unprotected, metal 
will corrode rapidly until failure occurs or until 
polarizing effects or shifting stresses permit reforming 
of the protective films. Corrosion will then by 
general until stress concentrations again rupture th 
protective film. 

In experimental confirmation of this theory, tensi 
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ns of fine-grained cartridge brass were ex- 


specin 
nosed under a tensile stress of approximately 45 
perce!!! of the vield strength in an ammonia, carbon 
dioxide, wate r vapor, and air atmosphere for periods 
ranging from 2 to 8 hr. Specime ns exposed for 2 to 


4 hr contained pits with their apexes at grain bound- 
aries, aS Is shown in figure 5, a and b. The first 
evidence of stress-corrosion cracks was found in 
specimens e xposed for 6 hr. These cracks extended 
from the roots of pits into the alloy along the grain 
boundaries in some of these specimens (fig. 5, ¢ and 
d). 

The residual stresses in the large-grained speci- 
mens of alpha brass were believed to be extremely 
small. The resolved tensile stresses on the planes of 
the stress-corrosion fractures were computed as 
fractions of the applied tensile stress and are given in 
table | It is seen that for most specimens the re- 
solved tensile stress was 85 percent or more of the 
applied stress. ‘Thus, the stress patterns were gen- 
erally in such directions as to produce stress concen- 
trations at the roots of notches formed by the 
corrosive process. ‘This satisfied one of the condi- 
tions for the propagation of stress-corrosion cracks. 

There was no appreciable difference between the 
results obtained on the cast and fully annealed 
wrought specimens of the 90-percent Cu-10-percent 
Zn alloy from the same heat. 


3.2. Beta Brass 


Eighteen crystals of beta brass (from four speci- 
mens) containing 47 transerystalline stress-corrosion 
cracks were studied. The normals to the crystallo- 
graphic planes on which these cracks were believed 
to have originated were determined and are plotted 
in a single stereographic triangle in figure 6. It is 
seen from these data that there is no plane or family 
of planes in beta brass that is particularly susceptible 
to stress-corrosion cracking. 

If the stress-corrosion cracking of beta 
associated with the rupture of the protective film, 
as is believed to be the case for alpha brass, it must 
be associated with either elastic or plastic deforma- 
tion of the material. Webb [78], Rinehart [79] and 
Good [20] have shown that single crystals of beta 
are very anisotropic in the elastic range. 
Young’s modulus is 8.2 times as great in the [111] 
direction as in the [100] direction, and the shear 
modulus is 12.8 times as great in the [100] direction 
as in the [111] direction. It was also shown that 
there is 2 linear relation bet ween | Kor LG and f(l,m, 
n), where f(/,m,n) = (Pm? *4t. 2/7) and/,m.and n 
are the direction cosines of the normals to any 
crystallographic plane with respect to the cubic 
axes. The values of 1/# and 1/@G have been com- 
puted for both the directions parallel to the stress 
axis and those norma! to the transerystalline stress- 
corrosion cracks studied and are given in table 
It is recognized that, because of the restraining 
effects of the neighboring crystals, the strain pro- 
duced in a given crystallographic direction in a 
crystal in a polycrystalline aggregate may be different 
from that calculated for a single crystal. 
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stress-corrosion cracks in beta brass. 


Ficure 6. Normals to 


Data obtained at the Bureau indicated, however, 
that there was greater elastic tensile strain for a 
given stress in polycrystalline grains of cast beta 
brass, oriented so that the [100] directions were closer 
to the tensile axis, than for grains with the [100] 
directions oriented 10° or more from that axis. 
Most of the erystals studied in this investigation, 
however, were so oriented that the [100] directions 
made angles greater than 10° with both the tensile 
axis and the normals to the planes of fracture. 
Hence, the possibility of microscopic plastic deforma- 
tion by slip was considered. The critical shear 
stress for slip may be defined as r=o cos ¢-cos X, 
where o is the applied stress, and @ and X are the 
angles between the direction of the applied stress 
and the slip plane and the slip direction (in the slip 
plane), respectively. Slip in beta brass has been 
reported [2/7] to occur on the (110) planes in the 
[111] direction. The angles that both the tensile 
axes and the normals to the planes, on which the 
stress-corrosion cracks were believed to originate, 
made with the various (110) planes and the [111] 
directions (lying in the respective (110) planes) were 
measured for each crack studied. The products of 
the cosines of these angles were determined for each 
pair of angles for each stress-corrosion crack, and 
the maximum value or values are given in table 2. 
The values obtained for the normals to the planes 
upon which stress-corrosion cracking originated are 
plotted in figure 7. There were also available orien- 
tation data for 81 large crystals in two polycrystalline 
beta brass tensile specimens. The values of the 
product COS @-COS X (with respect to the tensile axes 
of the specimens) were determined and are also 
plotted in figure 7. If the probability of slip in a 
given crystal incre ases as the value of cos @-cos A 
approaches 0.5, it isseen that there is a greater prob- 
ability of lie in the crystals containing stress-cor- 
rosion cracks than in the randomly selec ted crystals 
in the tensile specimens. 

A study of the data in table 2 indicates that in 
many cases both the probability of slip (cos ¢-cos ) 
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Values from table 2 are weighted so that they can be plotted on the same scak 
as data from randomly distributed crystals. @=angle between normal to slip 


plane and stress-corrosion crack; \=angle between the normal to stress-corrosion 
crack and the slip direction lying in the slip plane Critical resolved shear 
stress =o COS: @COSA ), cryst i] containing stress-corrosion cracks; X, randomly 


distributed crystals 


and the reciprocal of the shear modulus (shear strain 
per unit shear stress) are large not only for stresses 
parallel to the stress axis but also for those normal 
to the stress-corrosion crack. Generally, in other 
cases, if one of these values is small the other is 
large.* 

Smith [22] found experimentally that the grain- 
boundary energy of beta grains in an alpha-beta 
brass is only about three-quarters that of the alpha 


boundaries. If the grain-boundary energy in this 


alloy increases with increased difference in orienta- 
tion of the crystals, as has been shown for certain 
other metals, it may be expected that pitting will 


occur at certain high-energy boundaries more readily 
than on grain surfaces. Hence, some 
sion cracks may be expected to originate at grain 
boundaries in the beta brass by the same process 
described for alpha brass. Metallographic examina- 
tions at the Bureau showed that a small percentage 
of the stress-corrosion cracks do originate at grain 
boundaries. However, these become transcrystal- 
line before they have penetrated to the depth of a 


stress-corro- 


full grain into the metal. Evidently factors other 
than high-energy grain boundaries determine the 
4In erystal 9 both cos ¢ cos A and 1/G are small; however, 1/E has a value of 


dyneX 10 It is suggested that in this case the film rupture resulted 


Strain 
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from elastic tensil 








paths followed by stress-corrosion cracks in beta 
brass. 
It is postulated that regions of stress concen! ration 


produced by pitting corrosion (or in a few instances 
by grain- -boundary attack) in beta brass nucleate 
microscopic slip in the suitably oriented ¢1 stals. 
This slip ruptures the protective film (formed in the 
atmosphere or in the corroding medium) exposing 
narrow areas of film-free metal that are anodic to 
the film-covered metal, and_ stress-corrosion then 
progresses as was indicated for alpha brass.  Iy 
further support of this thesis it should be noted 
that plastic deformation has been reported |23] and 
shown [24] to be associated with the transerystalline 
stress-corrosion cracking of magnesium alloys. Ip 
the Bureau, slip lines have recently been found 
associated with transcrystalline — stress-corrosion 
cracking of type 304 stainless steel. 


4. Summary 


The crystallographic planes followed by stress. 
corrosion cracks in beta brass and the angle between 
the planes in the crystals on either side of the inter- 
crystalline cracks in alpha brass were determined by 
» ray diffraction methods for specimens of beta 
brass and two compositions of alpha brass. 

The crystallographic planes bounding intererys- 
talline stress-corrosion cracks in alpha brass and on 
which the stress-corrosion cracks originated in beta 
brass were randomly distributed. There was no 
plane or family of planes that was particularly 
susceptible to stress-corrosion cracking. 

About 70 percent of the stress-corrosion cracks 
studied in alpha brass originated between crystals 
in which the bounding planes were oriented at angles 
of 13° to 23° to each other. This is the range in 
which the maximum interfacial energies have been 
reported for other materials. It is suggested that in 
a corroding medium the grain boundaries having the 
highest energies become the loci of pits. In the 
absence of stress these pits are most probably self 
limiting, but with the application of favorably 
oriented tensile stresses they act as areas of stress 
concentration that may produce rupture of the 
normal protective film at the roots of the pits. This 
film-free metal is rapidly attacked and early inter- 
crystalline failure of the metal may result. 

Assuming that there is a critical resolved shea 
stress for slip in beta brass, it was shown that there 
was a greater probability of slip occurring in crystals 
containing transcrystalline  stress-corrosion cracks 
than in randomly distributed crystals. It is postu- 
lated that in a corrosive medium, microscopic slip, 
most probably nucleated by pitting, ruptures the 
protective film on the surface of the metal. The 
film-free metal is rapidly attacked and transerystal- 
line failure of the brass results. 

The resolved tensile stresses normal to the stress- 
corrosion fractures were generally 80 percent or more 
of the applied tensile stresses in both the alpha and 
beta brasses, thus indicating that these stresses were 
favorably oriented for the propagation of 
corrosion cracks. 


stress- 
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bets TABLE 2. Probability of slip, f(cos @, cos X) and reciprocals of elastic and shear moduli of stresses applied along the stress axis 
a and normal to transcr ystalline stress-corrosion cracks in beta brass. f(cos @, COS A) =CO8 @ COS \ where dg and Xarea ngles between 
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